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1.0 SUMMARY 


It is well known that there exists systematic differences in the fan stage noise between static 
(ground test) and flight conditions. These differences are significant in the current high 
bypass ratio engines with single stage, non-IGV fans. They have been designed to provide 
complete cutoff at the blade passing frequency (BPF) throughout the approach power 
range. However, results from static tests with these engines reveal a strong presence of the 
BPF tone. This tone has been associated with the steady and unsteady distortions in the 
static fan inflow field in particular the atmospheric turbulence, the ground vortex, test 
stand induced distortions and the interference of all these disturbances with the inlet 
boundary layer. Fan noise data from flight tests with these engines show significantly lower 
BPF levels, but usually indicate that complete elimination has not been achieved. The 
remaining fan tone has been associated with the atmospheric turbulence, inlet angle of 
attack related distortions as well as incomplete cutoff of rotor stator interaction noise. 

The primary purpose of the present investigation is 

• to define an atmospheric turbulence model 

• to identify the turbulence length scales that dominate in the generation of fan tone 
noise due to convected turbulence, and 

• to evaluate the difference in fan tone noise due to atmospheric turbulence between 
static and flight conditions. 

The investigation is based on analytical and empirical models and is divided into two segments. 
The first segment is concerned with the definition of the turbulence at the fan face, which is 
related to the atmospheric turbulence in the ambient flow field and is affected by the flow 
contraction between the ambient and the fan face conditions. The second segment is 
concerned with the fan noise generation due to convected turbulence. 

Based on a literature search an atmospheric turbulence model has been selected that predicts 
turbulence intensities, integral scales and spectra as a function of the altitude above the 
ground, the mean wind velocity, the surface roughness and atmospheric stability. An 
existing flow contraction model has been extended to provide all the one-dimensional spectra 
of the post-contraction turbulence for simple, isotropic pre-contraction turbulence. Narrow- 
band fan noise resulting from ingested anisotropic but homogeneous turbulence has been 
studied with existing fan noise models with particular emphasis on fan noise resulting from 
atmospheric turbulence in static and flight conditions. The primary results of the study are: 

• The difference in fan tone power at BPF due to ingested atmospheric turbulence between 
typical static test conditions and typical landing approach fliglit conditions is in the 
order of 30 dB. Based on this difference it is concluded that fan noise due to ingested 
atmospheric turbulence is negligible in landing approach conditions for the currently 
used high bypass ratio engines. 



Based on the analytical results it is concluded that fan noise due to atmospheric turbu- 
lence is below actual fan noise in flight conditions. It is therefore assumed that it is not 
necessary to simulate in-flight turbulence in static tests. However, fan noise due to 
atmospheric turbulence in static tests should be reduced to levels below broadband fan 
noise in flight conditions. An effective way to achieve this, is the reduction of the 
turbulence intensities in the static inflow field. The following guidelines are recom- 
mended for static test conditions: 

. usage of an inflow control device 

• low limits on wind velocities 

. small surface roughness scale in the test stand environment 

• large height of engine axis above ground 

. mean wind direction from the forward arc to eliminate reingestion 

Proper test stand and bellmouth inlet design and possibly inlet boundary layer suction 
are necessary to achieve this reduction. 

Rotor tone noise due to convected turbulence is dominated by a small range of transverse 
scales. For typical turbulence energy distributions this range covers about one decade 
and is centered around a transverse scale in the order of 30% of the rotor blade spacing 
at the blade tip. For a given turbulence variance, maximum sound power is achieved 
at transverse integral scales in the order of 25% of the rotor blade spacing at the blade 
tip. For a JT9D fan this is about 4. 1 cm and for the JTl 5D fan it is about 1 .3 cm. 

The difference in fan tone noise due to atmospheric turbulence between static and flight 
conditions is primarily due to the difference in the intensities of the turbulence involved 
in the fan tone noise generation. In static conditions the atmospheric turbulence is 
affected by a large inflow contraction, which results in: 

. large initial scale and therefore high initial turbulence energy level 
. large increase in the standard deviation of the transverse turbulence component 

• large streamwise integral scale at the fan face 

In the flight conditions the turbulent energy at large wavenumbers decreases with 
altitude, due to the nearly constant turbulence variance and the increase of the integral 
scales with altitude. 

The effect of inlet angle of attack on fan noise has been evaluated with data acquired 
during a fliglit test with a 747 aircraft powered with JT9D engines. Sound pressure 
measurements were made on the inlet wall, the fan duct wall and the fuselage surface 
over the full operational range of aircraft angle of attack. The measurements indicate 
that there is no obvious correlation between inlet angle of attack changes and changes 
in broadband noise as well as narrowband noise at BPF and its second and third harmonic. 
A circumferential variation of the buzzsaw sound pressure field in the inlet that increases 
with angle of attack, has been observed. 



2.0 INTRODUCTION 


By now it is well known that systematic differences exist between the fan stage noise in 
static and flight conditions (Ref 1 to 6). These differences became apparent with the intro- 
duction of the current high bypass, single stage, non-IGV fans. These engines are designed 
to provide complete cutoff at the fan blade passing frequency throughout the approach 
power range. However, results from static tests with these engines revealed a strong presence 
of the blade passing frequency (BPF) tone. Most of the corresponding flight test data show 
significantly lower BPF levels, but indicate that complete cutoff has not been achieved. The 
appearance of the BPF tone in static tests has been related to rotor noise due to steady and 
unsteady distortions in the fan inflow field, such as atmospheric turbulence, the ground 
vortex, test stand induced distortions and the interference of all these distortions with the 
inlet boundary layer (Ref 5 to 8). Much progress has been made in the simulation of flight 
inflow conditions with inflow control devices (Ref 4, 9 and 10). Further improvements have 
been achieved with the combined application of inflow control and inlet boundary layer 
suction (Ref 1 1 ). Another approach that has been used successfully to simulate in-flight 
conditions is the evaluation of scaled fan models in wind tunnels (Refs 10, 12 and 13). 

The cause of the remaining fan tone in flight is not known. It has been associated with 
atmospheric turbulence, inlet angle of attack induced distortions as well as incomplete 
cutoff of rotor-stator (primary or secondary) interaction noise. Little work has been done 
in the evaluation of fan noise due to atmosplieric turbulence in flight conditions. A Boeing 
in-house report used simple turbulence and fan noise models and showed significant 
differences between the two conditions. Substantial differences in the unsteady pressure 
distributions on propellor blades between static and flight conditions, which have been 
related to convected atmospheric turbulence, are reported in Reference 15. 

It is the primary purpose of the present study to evaluate the difference in fan noise due 
to atmospheric turbulence between static and flight approach conditions. The investiga- 
tion is divided into two segments. Tlie first segment is concerned with the definition of 
the turbulence at the fan face, which is related to the atmospheric turbulence in the ambient 
flow field and is affected by the flow contraction between the ambient and the fan face 
condition. An existing model that describes the properties of the atmospheric turbulence 
(Ref 16) is used to define the turbulence in the inflow streamtube ingested by the fan. 

The effects of the inflow contraction between the ambient and the fan face conditions on 
the atmospheric turbulence are estimated based on an existing model (Ref 17), which has 
been extended to provide more detailed information. In order to compute fan noise resulting 
from convected turbulence with the selected fan noise models, it was necessary to represent 
the turbulence at the fan face by a distribution of discrete distortion elements. A model 
that computes the statistical averages (spectra, correlations, length scales and variances) of 
such distributions of distortion elements has been developed and is described. 

The second segment of the study is concerned with the fan noise generation due to convected 
turbulence. The range of transverse scales involved in fan noise generation is determined and 
a limited parametric study of fan noise due to convected turbulence is conducted. All the 
models mentioned previously are then combined to estimate the difference in fan narrowband 
noise due to atmospheric turbulence between static and flight conditions. Results based on 
measured and predicted turbulence at the fan face are compared. 


3 



An additional segment of the study is concerned with the inlet angle of attack effects on 
fan noise. Results from a flight test, conducted jointly by P&WA and Boeing with a 747 
aircraft powered by JT9D engines are presented. 

Use of trade names or names of manufacturers in this report does not constitute an official 
endorsement of such products or manufacturers, either expressed or implied, by the National 
Aeronautics and Space Administration. 

The author would like to acknowledge the technical support received throughout this project 
from T. L. Clark. 
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3.0 SYMBOLS AND ABBREVIATIONS 


a 

^mnk 

B 

c 

Ch 


C; 


fi 


G(k) 

h* 

u 


V 

^mn 

C 


Transverse scale of distortion elements 

Mode amplitude of spinning order m, radial order n and harmonic order k 
Half bandwidth used in the evaluation of narrowband tones 
Velocity scale of largest distortion elements 
Speed of sound 

Velocity scale of distortion elements 
Rotor blade chord 
Rotor blade section lift coefficient 
Specific heat of air at constant pressure 

Correlation coefficient of the i^^ turbulence velocity component 
Boundary layer thickness 

Coefficient, 1. for Gaussian eddy field, e/2 for vortex field 
Eddy or vortex volume 

tVi 

One dimensional spectrum of the i'-^ turbulence velocity component 
Gravitational acceleration 
Spectrum function 

Airfoil acoustic response function 

Modified Bessel function of the first kind and order C* 

Axial acoustic energy tlux 
Von Karman constant (.4) 

Wavenumber vector in turbulence velocity field 
Harmonic order in fan narrowband sound pressure field 
Axial wavenumber component of fan narrowband sound pressure field 
Mixing length 

Streamwise scale of distortion elements 


Contraction ratio 



Monin-Obukov length scale 
Circumferential distortion order 

Integral length scale of turbulence component along coordinate 
Exponent in /? probability density distribution 
Spinning order 

Axial mean flow Mach number 

Mean flow Mach number relative to fan rotor blade 

Mechanical fan tip Mach number 

Exponent in probability density distribution 

Radial order 

Constant in spectrum function 

Density of distortion elements 

Number of rotor blades 

Static pressure 

Acoustic pressure 

Probability density distribution 

Velocity amplitude of turbulence wave 

Occurrence rate 

BL/M 

Radial coordinate within distortion elements 

Radial location of distortion element center 

Radius of outer annulus wall 

Radius of inflow streamtube 

Correlation function of i^^ turbulence component 

Normalized duct eigen-function 

Normalized random variable 


Normalized random variable 



T 

u 


U 

V 

w 

W 

X 

z 


a 


a, (3 


7 

7 

r 

5 


e 

1 ? 

Jf 

X 

^mn 

P 


Static temperature 
Turbulence velocity component 
Streamwise turbulence velocity component 
Friction velocity 
Mean flow velocity 

Transverse turbulence velocity component 

Vertical turbulence-velocity component 

Acoustic sound power 

Cartesian coordinate 

Vertical coordinate 

Axial duct coordinate 

Surface roughness scale 

Height where turbulence variances become equal 

Rotor blade angle of attack 

Exponents in length-velocity scale relationship 

Constant in spectrum function 

Stagger angle 

Gamma function 

Three dimensional turbulence spectrum 
Delta function 

Contraction parameter ^£2 I ^ 

Hub to tip ratio 

Post-contraction wavenumber vector 
Turbulence wavelength 
Annular duct eigenvalue 
Separation vector 
Radial duct coordinate 



a 

a 

T 

0^0 

a 

subscripts 

o 

1 , 2,3 

a 

ij,k 

fi 

u,v,w 

A 

AP 

B 

E 

F 

P 

REF 

S 

SL 


Standard density of air 
Turbulence velocity standard deviation 
Density ratio 
Time delay 

Angular duct coordinate 

Angular eddy center coordinate 

Universal functions in atmospheric turbulence model 

Angular frequency 

Blade passing frequency 

Vorticity vector 


Ground surface condition 
Standard condition 

Indices for vector or matrix components 
Associated with transverse distortion scale 
Indices for vector or matrix components 
Associated with streamwise distortion scale 
Associated with turbulence velocity components 
Pre-contraction condition 
Airplane 

Post-contraction condition 
Gaussian eddy 
Approach flight condition 
Pulse related parameter 
Reference condition 
Static test condition 
Single streamline 
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ss 


Single scale 
Vortex 


V 

superscripts 

A Pre-contraction condition 

B Post-contraction condition 

->• Vector 

abbreviations 

erf Error function 

BMT Blade mounted transducers 

BPF Blade passing frequency 

PDD Probability density distribution 

Ri Richardson number 

Proportionality 

= Approximate equality 



4.0 TURBULENCE MODELS 


A major goal of this study is to provide a comparison of fan tone noise resulting from ingested 
atmospheric turbulence between typical static and flight test conditions (Sec. 5.3). This 
requires the definition of turbulence models that allow the description of the atmospheric 
turbulence at the fan face. A model that defines the statistical averages for atmospheric 
turbulence in typical test conditions is described in Section 4.1. The description of a model 
that estimates the effects of the flow contraction between the ambient and fan face condi- 
tions on the atmospheric turbulence is included in Section 4.2. In the evaluation of fan tone 
noise due to turbulence, the latter is represented by a homogeneous distribution of distortion 
elements. A model that allows such a representation of the atmospheric turbulence at the 
fan face is described in Section 4.3. 

4. 1 ATMOSPHERIC TURBULENCE MODEL 


4.1.1 BACKGROUND 

The planetary or atmospheric boundary layer extends between the Earth’s surface and the 
outer atmosphere. It absorbs the shear between the surface, with its no-slip condition and 
the geostrophic wind at the upper edge. It is characterized by large scales; its thickness is 
in the order of 600 meters. It normally is divided into 2 layers, the outer layer and the 
surface layer. In the latter, the flow is assumed 2 dimensional, unaffected by the coriolis 
force. The random motion within this boundary layer, atmospheric turbulence, is also 
characterized by its large scale. The separation between variations in the mean wind velocity 
and turbulence fluctuations is best described by the Van der Hoven spectrum, Figure 1 . 

Atmospheric turbulence for two conditions, static and flight approach power noise tests, 
is considered in the present study. In static tests, the atmospheric turbulence in the test 
stand environment, 0 to 20 meters above the ground is of concern. In flight tests, the 
turbulence at the approach path altitude (120 meters) above the noise certification location 
must be considered. Both of these conditions occur within the atmospheric boundary layer. 

Only a limited amount of data describing atmospheric turbulence at the fan face of statically 
operated engines is available. References 19 and 20. No measurements have been made in the 
corresponding turbulence field in flight conditions. Such measurements would be very 
difficult to make, due to the extreme noise environment in the engine inlet and the low 
turbulence intensities associated with the small scales that are of interest. This lack of 
experimental data makes it necessary to base this study on empirical and analytical models. 
Several models are available that provide the statistical averages describing atmospheric 
turbulence. References 21 througli 23, and extensive bibliographies are presented in 
References 23 and 24. The model selected for the present study is summarized in the 
following paragraph. To a large extent, it is based on information developed in Reference 25, 
and it is described in great detail in Reference 1 6. 
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4.1.2 MODEL DESCRIPTION 


Both test conditions in this study are limited to low wind velocities. The limits are imposed 
to assure high quality for the acoustic data and satisfactory engine operation in static tests. 

For the flight condition FAR Part 36 prescribes a 5.14 m/s (10 kn) upper limit on the mean 
wind velocity at 10 m above the ground. The limits for static test conditions are not 
regulated, but similar wind limits are used tliroughout the industry. They are normally 
specified at the altitude of the engine axis above the ground as a function of the angle 
between the wind direction and the engine axis. The limits are typically 2.68 m/s (6 MPH) 
in the forward arc and fall off to zero for crosswinds and winds from the aft arc. 

Turbulence in the atmospheric boundary layer is generated by two mechanisms, buoyant 
convection and shear stresses. The shear stress related turbulence production is proportional 
to the mean velocity gradient and therefore decreases with altitude, whereas the buoyant 
production is assumed invariant with altitude. For this reason the turbulence production 
due to shear stresses dominates in high winds and at low altitude. But in low wind velocities 
and at high altitude, the buoyant production can become significant. For this reason one 
of the requirements for the atmospheric turbulence model was the inclusion of the buoyant 
production terms, to properly simulate low wind velocity conditions. Another requirement 
was the proper description of turbulence at large wavenumbers, where the turbulence becomes 
isotropic. This is important in the evaluation of in-flight turbulence, where only a limited, 
high wavenumber segment of the total turbulence spectrum is of interest. Based on these 
requirements the atmospheric turbulence model described in Reference 16 has been selected 
for this study. 

In the lower layer of the atmospheric boundary layer (the surface layer), the mean flow field 
is assumed to be two dimensional. This layer covers the lowest 10 to 20% of the atmospheric 
boundary layer. It includes the static test environment and its upper boundary is close to 
the flight condition. It is also called the constant stress layer, due to the fact that the shear 
stress T remains approximately constant throughout it. 

T = p Ll.^“ (4. 1.01) 


with 

p = air density 

= friction velocity 


Outside of the viscous sublayer the shear stress is dominated by the Reynolds stress 


r = puw = p 



(4.1.02) 


with 

u, w = turbulence velocity components 
£ = mixing length 

U = mean flow velocity 


The coordinate system is shown in Figure 2. 
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Figure 1 .—Schematic Spectrum of Wind Speed Near the Ground 
From a Study of Van Der Hoven ( 1957) (From Ref 18) 


z 



Figure 2.— Coordinate System for Atmospheric Boundary Layer Mode! 
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If it is further assumed that the mixing length 2 is proportional to the coordinate z, then 
the above equations lead to the well-known logarithmic velocity profile 


with 

k 


z 


o 



J2 = 

0Z 



k . z 




— 


(4.1.03) 

kz 

"^^0 0Z 


z 


(4.1.04) 


von Karman constant = .4 
roughness scale 


The friction velocity u^ is the velocity scale in the surface layer and the standard deviations 
of the turbulence components are proportional to it. Equation (4.1.04) is used to determine 
based on the wind velocity at a reference altitude zj^£p and the roughness scale z^. 


~ U (zref) 


k — 
Cn 


1 


(4.1.05) 


The friction velocity and therefore also the turbulence standard deviations are proportional 
to the reference wind velocity and increase with increasing surface roughness scale. The 
latter is a measure of the surface roughness and is in the order of 1 0% of the size of the 
roughness elements. Representative roughness scales are shown in Figure 3. In Reference 16, 
a roughness scale of .05 meters has been suggested for the airport environment. 

Equation (4. 1 .04) is based on the assumption of constant shear stress. It loses its validity 
when z approaches the size of the roughness elements and is therefore limited to conditions 
where z » Zq. 

Thermal Stability 


So far only turbulence produced by shear stress has been considered. In a neutral atmosphere 
the pressure p and temperature T follow the isentropic relationship and buoyant produc- 
tion is zero. 3 p 


• ■ 

9z Cp • p 9z Cp 

with 

g = gravitational acceleration 

Cp = specific heat at constant pressure 


(4.1.06) 
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Terrain description of area within several kilometers upwind of site 
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[■ Centers of cities with very tall buildings 


Centers of large towns, cities 

Centers of small towns 
Outskirts of towns 


[- Many trees, hedges, few buildings 


Very hilly or mountainous areas 


Forests 


Fairly level wooded country 




Many hedges 

Few trees, summer time 


Isolated trees 
Uncut grass 


Few trees, winter time 
Cut grass (^3 cm) 


Natural snow surface (farmland) 


Off-sea wind in coastal areas 


|- Farmland Long grass (^60 cm), crops 
Airports (runway area) 

I- Fairly level grass plains 


/ 




Calm open sea 


Desert (flat) 
(- Large expanses of water 


Snow-covered flat or 
rolling ground 


Ice, mud flats 


Typical 

static 

test 

environment 


Typical 
L flight 
test 

environment 




Figure 3. — Values for Surface Roughness Scale Zq (From Ref 22) 
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If however, the temperature gradient deviates from the isentropic relationship, buoyant 
production can either increase or decrease the stress produced turbulence. The atmospheric 
stability and its effect on the atmospheric turbulence are described by stability parameters. 
The best known parameters are the Richardson number (Ri) and the Monin-Obukov length 
scale (£'). (Ri) and (z/S') represent the production ratio between buoyant and stress generated 
turbulence. The parameter used in the selected model, the length scale £' is described in 
Reference 25. 


£' = 



(4.1.07) 


The sum in the denominator represents the difference between the actual and the isentropic 
temperature lapse rates. £' can be determined with simple mean velocity and temperature 
measurements. The following conditions can be defined based on the actual temperature 
lapse rate. 


unstable atmosphere 




.0 


neutral atmosphere 


n =_A 

3z Cp 



.0 


stable atmosphere 


^ >_A 

3z Cp 



.0 


Mean Wind Velocity Profile 


Turbulence generated by buoyant convection affects the mean velocity profile. For non- 
neutral conditions equation (4.1.03) that defines the wind velocity gradient becomes 


oz 



(4.1.08) 


with (z/£') a universal function. 
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So far the shear stress has been assumed constant, independent of altitude. But, in reality, 
it decreases with altitude and disappears at the edge of the boundary layer. In the selected 
model the shear stress is represented by the following relationship 


T 


= P = -To 




(4.1.09) 


with 

Tq = shear stress at surface 
u.^ = friction velocity at surface 

d = boundary layer thickness 


Equation (4.1.08) then becomes 


U = 




(4.1.10) 


(4.1.11) 


With <f >2 and 03 universal functions of (z/i2') which are shown in Figures 4 and 5. 

The surface friction velocity is determined from the mean wind velocity at the 
reference altitude the length scale 9! and the surface roughness scale z^. 

k . u(^REf) 

^REF'^M /^REF\ ^ref 



^3 


(7^) 


(4.1.12) 


Turbulence Standard Deviation 

The most consistent information concerning standard deviations of atmospheric turbulence 
is available for the vertical component. This is partially due to the lack of large scale motion 
close to the ground for this component. According to similarity analysis, the standard 
deviations of atmospheric turbulence are proportional to the friction velocity u^ and a 
function of the stability parameter (z/C'). 
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(4.1.13) 


The universal function 04 is shown in Figure 6. It assumes a value of 1 .3 for neutral 
conditions and increases with decreasing stability. For highly unstable conditions, it 
approaches a relationship characteristic for a free convection layer. For stable conditions, 
it decreases slowly, and above the critical value of 1., it drops to zero indicating that 
turbulence cannot exist above this stability level. 

According to equations (4.1.12) and (4.1.13) the standard deviation of the vertical turbulence 
component is proportional to the reference velocity and increases with the surface roughness 
scale and decreasing stability. For neutral stability it is insensitive to altitude. 

Close to the ground, the standard deviations of the horizontal turbulence components deviate 
significantly from the one of the vertical component. This is mainly due to the fact that 
large scale motion for those components is not suppressed as much close to the ground. 

In the selected model, it has been assumed that the standard deviations for the two horizontal 
turbulence components are the same. They are proportional to and the proportionality 
factor is a function of the ratio between z and the altitude zj, where the standard deviations 
of all 3 components become equivalent. 

z<zj (4.1.14) 


i ^ V 

\.177 + .823 —) 


o . , = u,, = (J. 


This empirical relationship is shown in Figure 7. Close to the ground the standard deviations 
of the horizontal components are twice as large as the one of the vertical component. 

Integral Length Scales 

Most integral length scales are based on autocorrelations from single point measurements in 
the atmospheric boundary layer. Taylor’s hypothesis is used to compute integral length 
scales from integral temporal scales and the mean flow velocity. These measurements provide 
only information about the integral scales in the streamwise direction. As in the case of the 
standard deviation of the turbulence velocities, the data for the integral scale of the vertical 
component are the most consistent. Similarity analysis suggests that the integral length 
scales are proportional to altitude with a proportionality factor that is a function of the 
stability parameter (C'/z). Measured data do however indicate that the correlation between 
atmospheric stability and integral scale is insignificant. Most empirical relationships are 
based on straight proportionality between length scales and height. The proportionality 
factors are in the order of .35 to .5, a value of .5 is used in the selected model 
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Figure 6. —Normalized Standard Deviation of the Vertical 
Turbulence Component (From Ref 16} 



Figure 7. —Standard Deviations of the Horizontal Turbulence 
Components (From Ref 16) 


1.22 





= .5 z 


(4.1.15) 


Measured data for the other two length scales, and show significant scatter and it is 
difficult to represent them satisfactorily by simple empirical relationships. In the selected 
model, they are computed based on the standard deviations, the integral scale and the 
assumption of isotropic turbulence at large wavenumbers. For von Karman spectra the 
above assumption leads to the following relationship for the streamwise integral length 
scales. Figure 8. 


2L^ = 

' w' 

(4.1.16) 

.5 z 

for z < zj 


.5 zj 

for z > Zj 



Turbulence Spectra 

Atmospheric turbulence is neither homogeneous nor isotropic. Its statistical averages 
are a function of the vertical coordinate and depend on the orientation of the coordinate 
system. It is however assumed that the turbulence is homogeneous in horizontal planes 
and stationary at a fixed location. This makes it possible to define local, statistical averages. 
Even though the atmospheric turbulence cannot be considered isotropic, its one dimensional 
spectra exhibit characteristics typical for isotropic turbulence. At high frequencies the 
spectra follow the inertial subrange relationship characterized by its (-5/3) exponent. In 
this frequency range, the ratio between the spectrum levels of the transverse and the stream- 
wise turbulence components becomes 4/3 (Ref 26), and the cross spectrum between the 
streamwise and the vertical turbulence components tends to zero (Ref 27). Based on these 
observations, it has been suggested that atmospheric turbulence is isotropic for wavelengths 
smaller than z/3 (Ref 21). Even though turbulence is not isotropic for wavelengths above 
this limit, it still is well represented by spectra developed for isotropic turbulence as long 
as the integral scales and variances are permitted to be different for the 3 turbulence compo- 
nents. The von Karman spectra satisfy the requirements for isotropic turbulence at large 
wavenumbers and fit experimental data well, especially in the inertial subrange. They are 
used to represent the spectral distribution of the turbulent energy in the selected turbulence 
model. 
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(4.1.18) 
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(4.1.19) 


The corresponding normalized spectra are shown in Figure 9. 

4.1.3 RESULTS 

The atmospheric turbulence model is used in Section 5.3.3 in the comparison of rotor tone 

noise due to atmospheric turbulence between static and flight conditions. Some relevant 

results are presented in that section. 

The following principal results can be inferred directly from the equations of the selected 

atmospheric turbulence model. 

• The standard deviations of the turbulence velocity components are proportional to the 
friction velocity, which is proportional to the reference wind velocity and increases 
with the surface roughness scale. 

• In the neutral atmosphere the standard deviation of the vertical component is nearly 
independent of altitude. In unstable conditions it increases, in stable conditions it 
decreases with altitude. 

• The integral scales of all the turbulence components increase with altitude. 

• Due to the fact that the standard deviation of the vertical turbulence component is 
nearly constant throughout the surface layer and because its integral length scale 
increases with altitude, its power spectrum level at large wave numbers decreases 
with increasing altitude. Figure 10. 
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Figure 10.— Spectra Based on the Same Turbulence Variance 
But Different Integral Scales L,- 
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4.2 INFLOW CONTRACTION MODEL 


4.2.1 BACKGROUND 

The most significant difference in the fan inflow field between static and flight conditions 
is the difference in the flow acceleration between the ambient and the fan face conditions. 

The fan face velocities are very similar for the two conditions, but the acceleration starts 
from the airplane forward velocity in flight, and from the local mean wind velocity in static 
tests. Similar flow contractions are present in wind tunnels between the settling chamber 
and the test section. It is primarily due to those, that aerodynamicists have been interested 
in the fluid mechanic characteristics of such contractions. It has been found that they have 
a significant effect on the freestream turbulence carried by the mean flow. This phenomenon 
has been studied with analytical models (Ref 28 to 30) as well as experimental investigations 
(Ref 31). Prandtl (Ref 28) studied the contraction effects on- turbulence with simple 
disturbance elements. In the evaluation of the streamwise disturbance velocity, he compared 
stream filaments that entered the flow contraction at different velocity levels. Assuming 
that the gain in kinetic energy due to the flow contraction is the same for both filaments, 
he concluded that the disturbance velocity decreased proportional to the inverse of the 
contraction ratio, which is the ratio between the post- and pre-contraction mean flow 
velocities. For the disturbance velocity normal to the mean flow he investigated a cylindrical 
element with the axis of rotation parallel to the mean flow. Based on the assumption of 
conservation of angular momentum, he found that the disturbance velocity increased propor- 
tional to the square root of the contraction ratio. 

A more sophisticated model was introduced by Taylor (Ref 29). He related pre- and post- 
contraction vorticity througli Cauchy’s transport equations for vorticity and introduced one 
of the key assumptions used in all contraction models: the relative motion between neighbor- 
ing particles due to turbulence must be small in comparison with the relative motion due to 
the flow contraction. Taylor found explicit solutions for simple, cellular vorticity distributions. 

Ribner and Tucker (Ref. 17) and Batchelor and Proudman (Ref. 30) used Taylor’s approach 
and investigated the effect of a flow contraction on isotropic turbulence. They showed that . 
for this type of turbulence the contraction effect on the turbulence variances is independent 
of the initial spectral distribution. Ribner and Tucker also investigated the effect of an 
axisymmetric contraction on initially isotropic turbulence with a spectral distribution based 
on the Dryden spectra. 

The turbulence model used in the present investigation is basically the one described in 
Reference 1 7. However, additional one-dimensional, post-contraction spectra have been 
derived. 


24 



4.2.2 REVIEW OF ANALYSIS 


An accelerating flow field stretches fluid elements in the streamwise direction and causes 
them to contract laterally. The sides of a cubical element change from D to C|D, 62^ and 
1230 as it travels from the pre-contraction position A to the post-contraction position B 
(Fig. 11). ^2 is defined as the mean velocity or contraction ratio and ^2 and £3 are the lateral 
contraction ratios. They are related through the continuity equation 

Ug Pg 

o> 2i • 1 £1 = a = — (4.2.01) 

^A ^A 


with ff the density ratio. 

A separation vector between two fluid particles will change from to Xg as it is convected 
from A to B, with 


^iB “ ®i • ^iA 


(4.2.02) 


Fluid filaments are therefore turned into the mean flow direction, unless they are aligned 
perfectly normal to the mean flow. 


The analysis described by Ribner and Tucker in Reference 17, is based on the following 
assumptions. 

• The relative motion between neighboring particles due to turbulence is negligible in 
comparison with the relative motion due to the flow contraction. 

• The contraction is rapid enough so that inertial effects dominate and viscous effects 
become negligible. 

• The lateral dimension of the flow field is large in comparison with the turbulence scale. 

• The mean flow field is one-dimensional. 


• Density fluctuations due to turbulence are negligible. 

Some of these assumptions are in conflict with each other. A rapid distortion for instance 
will provide one-dimensional flow only on the axis of the flow contraction. The require- 
ment of turbulence scales small in comparison with the lateral dimension of the flow con- 
traction will result in small scales which in turn may lead to turbulence dissipation in the 
contraction. But, despite these shortcomings, this analysis is one of the most appropriate 
tools to evaluate the inflow contraction effects on atmospheric turbulence. 

The analysis is based on the assumption of conservation of circulation. Cauchy’s transport 
equations for vorticity in Lagrangian form are used to define the post-contraction vorticity 
J2g based on the initial vorticity 12^ and the mean flow contraction. 
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Due to the assumption of one-dimensional mean flow and negligible relative motion due to 
turbulence, the distortion tensor reduces to its diagonal elements and the vorticity equations 
become 


9x-g 9x^ 

12;® = a 12;^ with 

' ' 9x 


iB 
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iA 


(4.2.04) 


These equations, together with the continuity equation, describe the contraction effect on 
turbulence. Due to the linearizing assumptions each turbulence wave (each Fourier component 
of the turbulence field) has to satisfy this system of equations. A Fourier component of the 
turbulence field at the pre-contraction location is represented by 

Ui^ = e^^^^ (4.2.05) 


The relative motion due to turbulence is negligible throughout the contraction and this 
turbulence wave is transformed into 




(4.2.06) 


at the exit of the contraction. The wavenumber vectors are related by 
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(4.2.07) 


See also Figure 1 1 . 


The amplitude vector Q for the Fourier components is determined with the continuity 
and vorticity equations 
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(4.2.08) 
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Based on these relationships Ribner and Tucker derived the equation for the components of 
the three-dimensional spectrum tensor For an axisymmetric contraction { 9.2 = 63), it 

assumes the following form: 
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( 4 . 2 . 09 ) 


( 4 . 2 . 10 ) 


The components of the three-dimensional spectrum tensor r^j^( 5 f) are expressed as a 
function of the contraction ratios and the corresponding pre-contraction tensor. Further 
simplifications are possible if the pre-contraction turbulence is isotropic. The pre- 
contraction spectra then assume the following form: 

Py A (k) = G(k) dy - kj kj ) (4.2.11) 

with 

9 jj =1. for i = j 
.0 fori^^j 

The spectrum function G (k), characterizes the spectral distribution of the turbulence. 

For such pre-contraction turbulence, the diagonal terms in the post-contraction spectrum 
tensor are 
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( 4 . 2 . 12 ) 


( 4 . 2 . 13 ) 
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F 33 ^ (?f) assumes the same form as F 22 (?f) since the initial turbulence is isotropic and the 
contraction is axisymmetric. 

Three-dimensional spectra provide extremely detailed information about the turbulence 
field, since the independent variable is a vector. It is, however, difficult to present them and 
takes an infinite number of standard two-point measurements to determine them experimentally. 

Standard one or two-point turbulence measurements provide only one-dimensional spectra, 
with a single wave-number component as the independent variable. They represent the power 
of all waves with a given wavenumber component independent of the magnitude of the other 
two wave number components. They are determined by integration of the three-dimensional 
spectra over two wavenumber components. 


OO 

(‘‘j) " ^ ff 

— OO 


(4.2.14) 


For a given spectrum function G (k), the one-dimensional post-contraction spectra can be 
determined by integration of the three-dimensional spectra given in equations (4.2.12) and 
(4.2.13). This integration is however only possible for simple spectrum functions. 

One dimensional spectra, that are frequently used to describe turbulence at large Reynolds 
numbers, are based on the following spectrum function 


G 



(4.2.15) 


where N and 7 are constants related to the turbulence variance and integral scale (Ref. 17). 



£_ 

L 



(4.2.16) 


with 

a = turbulence variance 
L = turbulence integral scale 


The one-dimensional spectra for the isotropic, pre-contraction turbulence based on this 
spectrum function are frequently referred to as the Dryden spectra. 
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(4.2.17) 
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(4.2.18) 
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The corresponding post-contraction spectra can be derived with equations (4.2.12) through 
(4.2.15). They form a matrix of nine spectra 


F]B (Ki) 

i^i) 

F 3 B (jfj) 


F,B (kj) 

F2® N 

F3® {^2} 


Fl® W 
F2® {^3) 
F 3 ® N 


It can be shown that all the corresponding cross spectra are zero. Due to the isotropic initial 
conditions and the axisymmetric contraction the following equivalences exist: 


Fl® N =F,B (JC3) 



F2 


B 





The entire matrix is therefore determined by five spectra. Two of them F} and 

F 2 ( 5 ^ 1 ) have been derived by Ribner and Tucker. They only provide information about 
the streamwise correlation of the turbulence. In the generation of rotor tone noise due to 
convected turbulence, the streamwise as well as the transverse correlations in the turbulence 
field are significant. For this reason, the remaining three spectra have been derived to 
complete the matrix. The derivation and the resulting equations are quite lengthy and are 
presented in the Appendix. 


The integral length scales at the post-contraction location are computed with the value of 
the corresponding spectrum at the wavenumber origin. 
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(4.2.19) 


The effect of an axisymmetric contraction on the variance of the initially isotropic turbulence 
is determined by integration of the three-dimensional spectra in equations (4.2. 12) and 
(4.2. 13) over all three wavenumber components. The results as derived by Ribner and 
Tucker are: 
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(4.2.20) 


(4.2.21) 


The correlation functions are determined by Fourier transformation of the one-dimensional 
spectra. 


OO 

Rj ^Xjj = J Fj cos kj Xj Skj 


(4.2.22) 


Correlation functions, as well as spectra, have been normalized with the turbulence 
variances and integral scales 
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(4.2.23) 
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(4.2.24) 


Equations (4.2.17) through (4.2.24) have been programmed and results are presented in 
the following section. The remaining integrals in the equations for the correlation functions 
are solved by numerical methods. 

4.2.3 RESULTS 

A limited parametric study has been conducted to evaluate the effects of the fan inflow field 
contraction on convected turbulence. The mean flow Mach numbers in this inflow field are 
reasonably small, (less than .6). The results indicate that for this range of Mach numbers the 
compressibility effects on the turbulence transformation are negligible. The results are 
therefore only a function of the contraction ratio, independent of the absolute Mach number 
levels. The results are presented as a function of the contraction ratio Cj, the mean flow 
velocity ratio. 
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Turbulence Standard Deviation 


The effects of the inflow contraction on the standard deviations of the turbulence components 
are shown in Figure 12. Results based on Prandtl’s and Ribner and Tucker’s equations are 
compared. They demonstrate the reasonably good agreement between the two models. 

Prandtl’s simple model seems to overestimate the effects of the inflow contraction. Both 
models show that the total turbulent energy increases with the contraction ratio. In the 
positive strain field of the flow contraction, the streamwise vorticity component absorbs 
energy from the mean flow field. The concurrent energy release from the other vorticity 
components, which are exposed to a negative strain field, is smaller and a net gain in turbulent 
energy results. Experimental and analytical data are compared in Figure 13. Uberoi’s data 
were measured in a wind tunnel flow contraction (Ref. 31). The engine inflow field data are 
based on turbulence measurements in the inflow field of a statically operated engine, outside 
of the inlet and upstream of the highlight plane (Ref. 32). The fan face data represent 
turbulence measurements made during a recent joint P&WA and Boeing noise reduction 
program. The turbulence was measured with split film anemometers at the fan face of a 
statically operated JT9D engine. The results demonstrate the increase in the separation between 
the standard deviations of the turbulence components with the contraction ratio. The 
discrepancy between the analytical and experimental results is mainly due to the insufficient 
attenuation of the streamwise turbulenee component. It has been suggested, that this is due 
to: 

• the tendency of turbulence to change towards isotropy 

• the decay of turbulent energy within the contraction 

• the contamination of the hot wire anemometer signal representing the streamwise 
turbulence component by acoustic signals and density fluctuations in highly anisotropic 
turbulence. 

Integral Length Scales 

The effect of a flow contraction on the integral scales is shown in Figure 1 4. The lines labeled 
Cj and £2 represent the streamwise and transverse contraction ratios and therefore represent 
the simple deformation model of Figure 1 1 . The integral scales are normalized with the scale 
parallel to the turbulence component in the initial isotropic condition. Depending on their 
orientation relative to the turbulence component, they assume a value of 1 . or .5 in the 
isotropic conditions. 

The integral scales L 21 and L 22 of the transverse turbulence component follow the corres- 
ponding contraction ratios quite closely. The third integral scale for the transverse component 
drops off very rapidly. This is due to significant negative values in the correlation function 
at large separation distances. 

The effects of the flow contraction on the integral scales of the streamwise component cannot 
be explained with the simple deformation model. The streamwise scale Lj j increases at a 
slower than expected rate and the transverse scale Lj2 increases rather than decreases with 
the contraction ratio. Both reflect the fact that the effects of the contraction on the turbulence 
are wave number dependent. 


32 


100 


10 





.01 h 


.001 U 



Ribner 

and 

Tucker 


= Streamwise contraction ratio 
a.j = Standard deviation of streamwise turbulence component 
O 2 = Standard deviation of transverse turbulence component 
0^= Standard deviation of precontraction turbulence 

=4 L_ 

1 10 100 


^1 


Figure 12.— Effect of Flow Contraction on Standard Deviation of Turbulence 
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Figure 13.— Ratio Between the Standard Deviations of the Turbuience 
Velocities in a Flow Contraction 
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One Dimensional Spectra 

The effect of the flow contraction on the one dimensional turbulence spectra is shown in 
Figures 15 through 18. The spectra in Figures 15 and 16 are normalized with the precon- 
traction variance and length scale and are presented as a function of the streamwise 
component of the precontraction wave number. The spectra of the transverse turbulence 
component in Figure 1 5 show an almost uniform energy increase (independent of wave 
number) with increasing contraction ratios. The energy increment seems to be wavenumber 
dependent only for contraction ratios smaller than about 3. The spectra for the streamwise 
component in Figure 1 6 demonstrate the significant decrease in the spectrum levels with 
increasing contraction ratio, as well as a significant change in the spectral distribution of the 
turbulent energy. This change in the spectral distribution is due to the energy transfer 
between the turbulence components. At large contraction ratios, it is apparent only in the 
spectra of the streamwise component, since its energy levels are much lower than the ones 
for the transverse components. 

Figure 17 shows the three one-dimensional post-contraction spectra for the transverse 
turbulence component. They are normalized with the post-contraction variances and 
integral length scales and show that the changes in the spectral distributions take place in the 
initial segment of the contraction (ratios less than 3.). The changes in the spectra F 2 (?f[) 
and F 2 (K 2 ) are minimal, but they are significant for F 2 Spectra such as the post- 

contraction spectrum F 2 (^( 3 ) are typical for unidirectional vortex fields as will be shown in 
Section 4.3. The corresponding spectra for the streamwise turbulence component in 
Figure 1 8 demonstrate that the spectral distribution changes up to the largest contraction 
ratios. These changes will however not be of much concern due to the fact that the streamwise 
component loses its energy and its significance decreases rapidly with increasing contraction 
ratios. 

Correlation Coefficients 

The three correlation coefficients for the transverse turbulence component are shown in 
Figure 19. The post-contraction coefficients represent results for contraction ratios larger 
than 3. The abscissa for the pre-contraction coefficients are all normalized with the same 
integral length scale measured parallel to the turbulence component. For the post- 
contraction coefficients it is normalized with integral length scales L 21 and L 22 - Only the third 
correlation coefficient is significantly affected by the contraction. It represents the transverse 
correlation in the direction normal to the turbulence velocity component and exhibits 
significant negative correlations at large separation distances. 

During a test program jointly conducted by P&WA and Boeing, turbulence data were acquired 
with split film anemometers in the inlet of a statically operated JT9D engine. Several 
sensors located on a radial strut measured the transverse circumferential turbulence component. 
Correlation between the different sensors provided information about the radial correlation 
of the turbulence. The results are presented in Figure 20 and exhibit large negative values at 
large radial separations. They therefore agree qualitatively with the predictions of the inflow 
contraction model. 
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Figure 16.— Effect of Flow Contraction on Spectra of Streamwise Turbulence Component 
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Figure 17.— Normalized Pre- and Post-Contraction Spectra for the Transverse Turbulence Component 
for Various Contraction Ratios C; 
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Figure 20.~Comparison of Measured and Predicted Post-Contraction 
Correlation Coefficients for the Transverse Turbulence 
Component 
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4.3 INFLOW TURBULENCE REPRESENTATION 


4.3.1 BACKGROUND 

The evaluation of blade passage frequency noise due to convected atmospheric turbulence is 
based on analytical models developed by T.L. Clark at The Boeing Company (Refs. 33 and 34). 
In these models, the convected turbulence is represented by a homogeneous distribution of 
distortion elements. The turbulence field at the fan face as predicted by the combined 
models of Sections 4. 1 and 4.2 is however presented in the form of one-dimensional spectra 
and correlations. To bridge this gap, an inflow turbulence model has been developed, that 
computes the one-dimensional spectra and correlations for random distributions of distortion 
elements. This model, which is described in this section, makes it possible to select the 
distributions of distortion elements that best simulate a given inflow turbulence field. 

The model is based on an approach similar to the one used by Melick in Reference 35. 

Melick’s model was developed for the study of compressor stall due to turbulence in the 
subsonic diffusor of supersonic inlets. The model attempts to simulate the turbulence 
resulting from the inlet boundary layer and shock-boundary layer interactions. It is based 
on a two-dimensional flow field and the distortion elements are axisymmetric vortices with 
their axis normal to the mean flow velocity. 

The turbulence model presented here describes the atmospheric turbulence at the fan face 
after it has been affected by a flow contraction. The distortion elements are axisymmetric 
with their axis parallel to the mean flow velocity. Their distribution is homogeneous in 
unbounded, three-dimensional space. The effects of inlet surfaces are neglected. 

4.3.2 ANALYSIS 

The acoustic models used in this study compute the tone noise resulting from distortion 
elements that are convected through a fan. In these models, it is assumed that only two 
unsteady velocity components contribute to fan tone noise, the streamwise component u 
and the transverse circumferential or angular component v. In the original report (Ref. 33) 
that describes these models, the distortion elements are represented as Gaussian eddies. The 
velocity distribution in such eddies is 

1 

u (x) = C| e 


(?) 


with 

C| = velocity scale 
£ = streamwise length scale 
a = transverse length scale 

(See also Fig. 2 1 ). These Gaussian eddies are used to describe unsteady velocities in the 
streamwise as well as the transverse direction. 
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In an attempt to provide a distortion element that is more realistic from the fluid dynamics 
point of view, a vortex element has been added to the acoustic models. It is used only to 
represent transverse velocity distortions, and its angular velocity distribution is 



(4.3.02) 


The velocity component U2 parallel to the X2 axis (Fig. 2 1 ) is therefore 


^ X3 ' 2 L a2 

V (x) = Cl e — e 

^ Q 


(4.3.03) 


The following analysis is based on the vortex element, but the results for the Gaussian eddy 
are very similar and are summarized at the end of this section. 

A vortex element convected along a streamline through point (0, X2, X3) will induce a velocity 
pulse at a point (0, ^2» ^3) 


X 3 - ^3 

with 

U = mean flow velocity 
t = time 

This velocity pulse is the same for all points in the (X2, X3) plane, except for its amplitude, 
which is a function of the relative position of streamline (x2, X3) and the point (^2> ^3)- ^ 
distribution of vortices convected along the streamline will induce velocity pulses at the 
origin (0, 0, 0) and at the point (0, ^2> ^3) (Figure 22). In Reference 36 it is shown that if 
the vortices are Poisson distributed along the streamline, the correlation function 
Ra (r, ^2» ^3) is equal to the product of the pulse correlation function Rp (t, ^2> ^ 3) and 
the average occurrence rate M. 
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(4.3.04) 


Ra ^ 2 ’ ^3) ~ M • Rp (r, ^3) (4.3.05) 


M = N . U . dx2 • dx3 


(4.3.06) 
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(4.3.07) 


with N = average number of eddies per volume unit. 


In a homogeneous distribution of vortices in unbounded space, Poisson distributed along all 
streamlines, the three-dimensional correlation function (t, ^ 3 ) can be computed by 

integration of (r, ^ 3 ) over X 2 and X 3 and a simple coordinate transform 


Rb 




(4.3.08 



(4.3.09) 


with V = variance. 

Rj^ is the correlation function for a distribution of vortices of equal scales a, C and Cj. It is 
assumed that in an actual turbulence field the random motion is due to a random distribution 
of vortices of various length and velocity scales. For this reason the scales a, 12 and c are 
assumed to be statistically independent random variables, each with its own one-dimensional 
probability density distribution. 


P(a), P(C), P(c) 

Beta probability density distributions are chosen for all three variables to provide flexibility 
(Figure 23); 


(4.3.10) 

with 

t = normalized independent variable (a/a^^gj^ ) 
r == gamma function. 
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To further enhance the flexibility of the model, the following deterministic relationship 
between velocity and length scales has been assumed. 






(4.3.11) 


In a random distribution of vortices, the density of vortices with a given set of scales is 

N . P(a) . P(C) . P(c) da dC dc (4.3. 1 2) 

The contribution of vortices of these scales to the correlation function is 

Rc(B = %(i) • P(a) . P(fi) . P(c) da d£ dc (4.3.1 3) 

The correlation function for a random vortex field is found by integration over the random 
variables. 
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(4.3.15) 
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(4.3.16) 


subscript V stands for vortex. 

The turbulence vari^ce is therefore proportional to the product of the mean square 
maximum velocity c2 of the largest vortices, the vortex density N and the volume E of 
the largest vortices. 



The corresponding one-dimensional correlation coefficients are found by setting 
ponents of the separation vector ($) to zero and normalizing with the variance v2 
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where and L2 are the corresponding integral length scales. 
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(4.3.21) 


m and n are the exponents in the Beta probability density distributions for the length 
scales a and £, while Sj^g^^ and upper bounds on the distributions of these 

variables. 

The normalized, one-sided, one-dimensional spectra are determined by Fourier transforma- 
tion of the corresponding correlation coefficients. 
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(4.2.24) 


Most of these equations are also valid for a random distribution of Gaussian eddies. But 
the velocity distribution in the Gaussian eddy is axisymmetric and therefore the one- 
dimensional correlation functions for the X 2 and X 3 components of the separation vector 
must be the same. 


The statistical averages for the random Gaussian eddy distribution are 
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4.3.3 RESULTS 


Equations (4.3.15) through (4.3.25) have been programmed. The integrals in the equations 
for the correlation coefficients and the spectra are solved by numerical methods. A limited 
parametric study has been conducted to determine the effects of the probability density 
distributions on the correlations and spectra and to find the appropriate distributions to 
represent the turbulence at the fan face. 

Turbulence Variance 

Equations (4.3.15) and 4.3.25) indicate that for given scale distributions P(a) and P(£) the 
turbulence variance is proportional to: 

N = density of eddies or vortices (number of eddies per unit volume) 

= mean squared maximum velocity of the largest eddies 
E = volume of the largest eddies. 

The volume E is of course related to the integral scales of the turbulence. The product 
Nc 2 appears only in the equation for the turbulence variance and the two variables cannot 
be separated. 

For the same set of parameters N, c and E and the same scale distributions the variance 
of a vortex field is larger by a factor of e/2 than the one of a Gaussian eddy field. 

Integral Length Scales 

The integral scales in a vortex field are the same as in a Gaussian eddy field, except for the 
transverse scale normal to the turbulence velocity component. The latter is equal to zero in 
a vortex field, whereas both transverse integral scales are the same in a Gaussian eddy field. 

The ratio between the integral and the maximum scale tends towards one for large values 
of m and towards zero for large values of n, with m and n the exponents in the probability 
density function. 

Spectra and Correlations 

Correlation coefficients and one-dimensional spectra for a field of identical vortices (a single 
scale vortex field) are shown in Figures 24 and 25. They do not compare well with the 
corresponding results presented in Section 4.2.3. The correlation coefficients are too wide 
at the apex and the spectra fall off too rapidly at large wavenumbers. Both discrepancies 
result from a lack of high wavenumber, small scale turbulence energy. 

Turbulence spectra that represent scale distributions rather than single scales are shown in 
Figures 26 through 31. They are based on equations (4.3.17) through (4.3.24). The 
parameters that control the shape of the spectra and correlations are the exponents mg, ng 
and /3 for the streamwise correlation and the corresponding spectrum and m^, n^ and a 
for the transverse correlations and the corresponding spectra. 
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The exponents rnj^ and /3, as well as and a appear only in the form of the differences 
m^-a and m^ - 13 in the equations for the spectra and correlations. The number of shape 
factors for each spectrum is therefore reduced to two. 

Spectra based on such scale distributions form an exponential relationship at large wave 
numbers. The exponent of this relationship is a function of the difference (mg - jS) or 
(m^ - a), see Figure 26. In actual turbulence spectra this exponential relationship represents 
the inertial subrange. In the von Karman spectra the exponent assumes a value of (-5/3), 
whereas in the spectra of Section 4.2 it is (-2). In order to simulate spectra with such small 
negative exponents, based on scale distributions with finite probability densities, it was 
necessary to introduce the exponents a and j3. The proper representation of the inertial 
subrange is important and the exponents (m^ - a.) and (mg - j3) are selected accordingly. 

The effect of the exponents n^ and ng on the spectra is shown in Figure 27. The variation 
of n is restricted by numerical limits, the values -.5 and 50. are close to the actual limits. 

The results indicate that even with the flexibility of the Beta probability density function, 
the spectra and correlations are strongly dominated by the corresponding functions of the 
single element. 

The normalized spectra and correlations for the vortex and the Gaussian eddy field are the 
same except for the ones in the direction normal to the mean flow and the turbulence velocity. 
For the Gaussian eddy field, the two transverse correlations and the corresponding spectra 
are the same due to its symmetry. 

Representation of the Turbulence at the Fan Face 

Based on the comparison of the spectra of this section with those of Section 4.2, it appears 
that in general the post-contraction turbulence field cannot be represented by simple 
Gaussian eddy and vortex distributions. However, in the evaluation of fan tone noise due 
to atmospheric turbulence, there are only two conditions of primary interest, static (ground 
test) and approach flight. The static condition is characterized by its large contraction ratios. 
The turbulence field at the fan face is dominated by the transverse turbulence component 
and can be well represented by a vortex field. The turbulence field associated with the 
streamwise component can be neglected. The approach flight condition is characterized by 
its small contraction ratios (between 1. and 2.). In Section 5.3 it is shown that in such 
situations the tone noise generation is dominated by the streamwise turbulence component. 

In this condition the turbulence at the fan face can be well represented by a distribution 
of Gaussian eddies. 

The results of Section 4.2 indicate that the exponential relationship at large wavenumbers 
is the same for the pre- and post-contraction spectra. In the large scale turbulence the 
exponent is (-5/3), representative for an inertial subrange. The same exponent has been 
selected for the representation of the turbulence at the fan face and the exponents m^, 
mg, a and j3 assume the following values 

7 4 

m3 = mj = .0 a = - (3 = j 
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Figure 27. —Normalized Spectra in Vortex Field for Various 

Exponents n in the Probability Density Distribution P(t) 







The exponents n^ and ng have only a minor effect on the spectra and are assumed to be 
zero 


na = ng = .0 

The normalized spectra for a Gaussian eddy field are compared with the von Karman spectra 
in Figure 28. They show reasonably good agreement. The normalized spectra of the vortex 
field based on the selected exponents are compared with the corresponding post-contraction 
spectra in Figure 29. Again the agreement is quite close. 

The contributions of various segments of the scale distribution to a spectrum is shown in 
Figure 30. This result demonstrates that in this particular representation of the turbulence 
the various segments of the spectrum are dominated by certain ranges of scales. The effect 
of the elimination of the largest scales within a distribution on the spectra is shown in 
Figure 31. For the assumed scale distribution (m = n = .0) the shape of the normalized 
spectrum does not change, only the integral scale and the turbulence variance are affected. 
This feature of the model is used in the evaluation of the atmospheric turbulence in flight 
conditions, where only a limited segment of the total range of scales is of interest. 
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Figure 28.— Comparison of Spectra Based on Von Karman Formulas 
and Gaussian Eddy Distribution 
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Figure 29.— Comparison of Spectra Based on Contraction Model 
and Vortex Distribution 
















5.0 ROTOR TONE NOISE STUDIES 


This section describes the analytical evaluation of fan tone noise resulting from atmospheric 
turbulence in the fan inflow field. The investigation is based on analytical fan noise models 
developed by T. L. Clark at The Boeing Company. The investigation is divided into 3 segments. 
In Section 5.1 tone noise resulting from single distortion elements is evaluated to determine 
the range of scales that cause significant tone noise levels. In Section 5.2 a parametric study 
of rotor tone noise due to convected turbulence is presented. The turbulence is based on 
predicted spectral distributions and is represented by a homogeneous distribution of distortion 
elements. In Section 5.3, the difference in rotor tone noise due to atmospheric turbulence 
between static and flight conditions is estimated based on the models and the results described 
in the preceding sections. 

5.1 ROTOR TONE NOISE DUE TO UNSTEADY DISTORTION 
5.1.1 REVIEW OF ANALYSIS 

This investigation is based on the analysis of Reference 33. The acoustic models described in 
that reference compute the sound pressure field at the harmonics of the blade passing 
frequency, radiated from a fan stage into an infinite, hardwalled annulus. The duct pressure 
field is described in the form of an eigenfunction expansion. 
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The analysis is based on the linearized theory of compressible fluids. Its applicability is 
therefore restricted to subsonic relative flow velocities. 


The noise generating mechanism considered in these models is the unsteady pressure force 
exerted by blade elements of the fan stage on the fluid. Such forces result from distortions 
in the mean flow field, which appear as unsteady velocity disturbances to the fan blades. 

The unsteady pressure distribution on the fan rotor blades is computed with the response 
function described by Naumann and Yeh (Ref. 37). Since this is a two-dimensional single 
airfoil response function, its validity is restricted to fan rotors of low solidity. The expansion 
coefficients, the mode amplitudes, are determined from the unsteady pressure distributions 
on the blades of the fan stage. 


The acoustic model used in the present study is the unsteady distortion or single eddy model. 

It computes the mode amplitudes resulting from a distortion element convected through a 
fan rotor. The considered elements are the same as in Section 4.3, the vortex and the Gaussian 
eddy. In the fan rotor plane, these distortion elements appear as distortions of 
modulated amplitude. The vortex distortion assumes the following form in the polar 
coordinates of the rotor plane. 
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with 

V = angular velocity component 
R, 4>, z .= coordinates of vortex center at t = 0 
a, £ = radial and axial vortex length scales 
cj = vortex velocity scale 

A similar equation describes the disturbance velocity field of the Gaussian eddy. 

After Fourier decomposition of its angular terms and Fourier transform of its time 
dependent term, equation (5.1.02) becomes 


(5.1.02) 
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(5.1.03) 


with 

C* = harmonic distortion order 

CO = distortion modulation frequency 

Ig* = modified Bessel function of first kind and order C* 

The harmonic of the disturbance velocity field is therefore composed of an infinite 
number of components, rotating at different angular velocities co/£* 


Vg*(p, 0, t) = F|£+(p) 


i(C*0-cot) 

J F2(co)e' 


(5.1.04) 


A fan rotor cutting through such a rotating distortion pattern creates a sound pressure field 
at frequencies separated by the modulation frequency from the blade passing frequency. 
The amplitude modulated velocity distortion in the plane of a rotating fan, representing a 
convected distortion element, therefore results in a sound pressure field in frequency bands 
centered around the blade passing frequencies. The modulation spectrum appears in the 
sidebands of these center frequencies. The equations for the mode amplitudes representing 
such a sound pressure field have been derived in Reference 33. 
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C* = m - kN 


(5.1.08) 


N = 
7? = 
Ch = 


8cl/3o; = 



7 - 
B = 


number of rotor blades 

hub to tip ratio 

rotor blade chord 

derivative of blade lift coefficient 

blade relative Mach number 

blade tip Mach number 

stagger angle 

half bandwidth 


In the derivation of the equation for the mode amplitude, it has been assumed that only 
the sound pressure field within a small frequency band is considered. This band is centered 
around the blade passing frequency or its harmonics and its width (2B) is assumed to be 
small in comparison with the blade passing frequency. It further has been assumed that all 
the modes of order m and n at frequencies within the band can be integrated into a single 
mode with amplitude I at the center frequency. The phase relationship between the 

modes that contribute to this narrowband mode varies with time and due to this the amplitude 
of the resulting mode varies with time. 

Equation (5.1.05) defines the mode amplitude resulting from a single vortex convected 
through a fan rotor. The corresponding equation for a Gaussian eddy is very similar and 
is shown in Reference 33. 
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The integrals in equation (5.1.05) are solved by numerical methods. The results are 
presented in the form of in-duct power levels W. The latter is determined by integration of 
the local, axial acoustic energy flux and is based on the following equation. 
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(5.1.09) 


with Rq the outer annulus radius. 

In the summation over the spinning and radial orders m and n the cutoff criteria is used to 
select the propagating modes. They satisfy the well known condition 
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(5.1.10) 


5.1.2 RESULTS 

A parametric study of rotor tone noise due to convected distortion elements has been 
conducted using the model described in the preceding section. The purpose of this study 
is to determine the transverse length scales of distortion elements, eddies or vortices, that 
cause significant tone noise levels. The rotor geometry represents a typical single stage, 
high bypass ratio fan (JT9D, CF6, RB.21 1) with 


hub to tip ratio 

r? = .38 

rotor solidity 

- 1.55 


The rotor blades are represented as twisted flat plates. The chord lines are aligned with the 
local mean flow velocity at the leading edge. 

relative mean flow angle (deg) tip 20 

hub 43 

The fan operating conditions cover the approach power range (high subsonic tip speed). 
The ranges of the parameters investigated in this study and the base case, the reference 
case are summarized in the following table. 



Parameter 

Range 

Base case 

Number of rotor blades 

12 46 

23 

Mechanical tip Mach number 

.7 .9 

.876 

Relative tip Mach number 

.745 .958 

.930 

Transverse eddy scale (duct radii) 

.01 1.0 

.1 

Streamwise eddy scale (duct radii) 

10. 

10. 

Radial location of eddy center (duct radii) 

.9 

.9 

Distortion velocity components 

streamwise. 

transverse 


transverse 

angular 

angular 

Distortion elements 

vortex, 

Gaussian eddy 

Gaussian eddy 


Radial locations of the distortion centers other than the 90% of the inlet radius have been 
evaluated during this study. But the results indicated that the sound pressure field due to 
distortions in the fan tip region is significantly higher than the one resulting from distortions 
at other radii. The results presented in Section 5.1 of this study are therefore limited to 
distortions at this particular radial location. The studies described in Sections 5.2 and 5.3 
are concerned with rotor tone noise due to convected homogeneous turbulence and are 
therefore concerned with distortions at all radial locations within the inlet. The parametric 
study was limited to the blade passing frequency and provided the following results. 

Transverse Distortion Scale 

The transverse scale controls the relative magnitude of the distortion velocity coefficients. 

A very large distortion with the scale in the order of the duct radius contributes primarily to 
the few lowest orders. For very small distortions with the scale in the order of one percent 
of the duct radius, the distortion coefficients are almost constant up to very high orders. 

In both cases, distortion coefficients of substantial magnitude contribute to non-propagating 
modes. It therefore can be assumed that there exists an intermediate range of scales where 
most of the significant distortion coefficients couple into propagating modes. The highest 
sound power levels are expected to occur in this intermediate range of scales. This is 
confirmed by the results presented in Figure 32, where the power level at the blade passing 
frequency is shown for the base case and a range of transverse scales. The results are based 
on the assumption that the maximum velocity is the same for all eddies independent of the 
transverse scale. For such eddies, the peak power level occurs at a transverse scale of about 
10% of the duct radius or 50% of the rotor blade spacing in the fan tip region. 

Distortion Velocity Components 

In most situations the unsteady blade lift forces are dominated by the effects of the distortion 
velocity component normal to the blade surface. Due to the small relative mean flow angle 
in the fan tip region, a streamwise distortion velocity component produces a larger blade 
normal component than the corresponding transverse distortion component. For the same 
velocity scale, Gaussian eddies representing streamwise distortions therefore produce higher 
sound power levels. For the investigated rotor geometry the difference in the fan tone power 
level between the two distortion velocity components is 7.6 dB for eddies in the fan tip region 
(Figure 32) and reduces to zero dB at the hub. Only the Gaussian eddy is used to represent 
transverse as well as streamwise distortion velocities. 
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Gaussian eddy 



Transverse scale in duct radii 


I \ I ! I I 

.05 .1 .2 .5 12 

Transverse scale in rotor blade tip spacings 


Figure 32.— Sound Power Level W as a Function of Transverse Scale 
for Various Distortion Elements 
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Distortion Elements, Gaussian Eddy and Vortex 

The distribution of the distortion velocity coefficients for the base case, a Gaussian eddy, 
and the corresponding vortex are shown in Figure 33. The relationship between the distortion 
coefficients and the harmonic order is very similar for the two elements at the dominating 
radial locations. However, the two elements differ significantly in the radial distribution of 
the distortion coefficients. For the Gaussian eddy, the radial distribution is a nearly even 
function compared to an almost pure uneven function for the vortex. Based on these 
distributions, it can be expected that for distortion elements in the fan tip region the eddies 
couple well into the lowest radial order modes, whereas the vortex couples better into the 
higher radial order modes. This is confirmed in Figure 34 which shows the modal power 
distribution for the two distortion elements. The total sound power is shown as a function 
of the transverse scale in Figure 32. Both distortion elements show similar trends, with the 
power levels falling off at very small as well as at very large transverse scales. The significantly 
lower levels for the vortex at intermediate transverse scales must be due to the fact that a 
substantial amount of the distortion is lost to higher radial order modes, which do not 
propagate. This comparison points out the significance of the radial velocity distribution, 
which is frequently neglected. 

Distortion Velocity Scale 

The assumption of a constant velocity scale independent of the transverse length scale does 
not conform with the energy distribution in actual turbulence. A simple relationship between 
the velocity scale and the transverse length scale can be used to approximate the turbulence 
energy relationship in the inertial subrange. It is assumed that each discrete eddy represents 
the eddies within a range of transverse scales, which is small in comparison with and propor- 
tional to the transverse eddy scale. Based on the results of Section 4.3 it can be shown that 
the turbulent energy AE within such a range is 


AE = Bi a2/3 

indicating that the energy is increasing with scale. In the same section it has been shown 
that the turbulence variance of a single scale eddy field is proportional to the product of 
the velocity scale squared and the eddy volume. In this simple representation, the eddy 
density and streamwise length scale are assumed to be independent of the transverse scale 
and therefore 


AE = B2 a^ 

These 2 equations for the turbulent energy provide the relationship between velocity scale 
and transverse length scale in this simple representation of turbulence by single eddies 

ci^=B3a“4/3 (5.1.11) 
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Figure 34.— Modal Power Distribution for Single Distortion Elements 



The assumption of a constant velocity scale corresponds to' an exponent of minus 3 in the 
inertial subrange, underestimating the energy in the small scales. The relationship of 
Equation (5.1.1 1) has been used to determine the relative sound power levels resulting from 
the eddies of different scales that represent the inertial subrange of a turbulence field. The 
following results are presented for the constant velocity scale as well as the relationship of 
Equation (5.1.1 1). 

Rotor Blade Gap in the Fan Tip Region 

This comparison is based on constant rotor solidity, the ratio between the total blading surface 
and the annular cross section. For a given rotor operating condition, the field of propagating 
modes at blade passing frequency is increasing with the number of rotor blades. Both the 
lower as well as the upper limit of the distortion orders contributing to propagating modes 
are increasing. It can therefore be expected that with increasing number of blades, the peak 
power level shifts towards smaller transverse scales. This is confirmed by the results shown 
in Figure 35. For the constant velocity scale data the peak occurs at transverse scales in 
the order of 50% of the blade gap in the fan tip region. For the case where the velocity scale 
is based on the turbulence energy distribution, it occurs at a transverse scale of about 25% 
of the blade gap. The following table summarizes the results for the JT9D and JTl 5D engines. 






Transverse scale at peak power level 




Number 







Inlet 

of rotor 


Constant 


Turbulence 


Engine 

radius 

blades 


velocity scale 

energy distribution 




a 

a N 

a 

a N 



RQ(cm) 

N 

— 

a(cm) 

— 


a(cm) 




Ro 

27T Rq 

Ro 

27tRo 


JT9D 

118. 

46 

.08 

.60 9.4 

.04 

.30 

4.7 

JT15D 

27. 

28 

.11 

.50 3.0 

.065 

.30 

1.8 


Mechanical Tip Mach Number 

With increasing tip Mach number, the field of propagating modes widens and distortion 
coefficients of lower, as well as higher harmonic orders start to contribute to propagating 
modes. High sound power levels will occur at larger as well as smaller transverse scales. 

This trend is demonstrated by the results shown in Figure 36. The results also indicate a 
shift of the peak power level towards larger transverse scales with increasing tip Mach number. 
The peak occurs at transverse scales in the order of 30% to 60% of the fan blade gap for the 
constant velocity scale and at 20% to 30% of the blade gap for the velocity scale based on 
the turbulence energy distribution representing an inertial subrange. 
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Rotor tip mach number = .874 

Gaussian eddy, transverse velocity component 




Figure 35.— Sound Power Level W as a Function of Transverse Scale 
for Various Numbers of Rotor Blades 
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23 rotor blades 

Gaussian eddy, transverse velocity component 



Transverse scale in rotor blade tip spacings 



.05 .1 .2 .5 12 

Transverse scale in rotor blade tip spacings 

Figure 36.— Sound Power Level W as a Function of Transverse Scale 
for Various Rotor Tip Mach Numbers 
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5.2 ROTOR TONE NOISE DUE TO TURBULENCE 


5.2.1 REVIEW OF ANALYSIS 

The turbulence rotor tone noise model is an extension of the unsteady distortion model 
described in the preceding section. It computes the mean power level at the blade passing 
frequency and its harmonics resulting from a homogeneous turbulence field convecting through 
a fan. The turbulence is represented by a homogeneous distribution of distortion elements 
(eddies) as described in Section 4.3. 

A convected eddy can be represented by an amplitude modulated distortion in the fan rotor 
plane. The sound pressure field resulting from a fan rotor cutting through such a distortion 
is not composed of pure tones, but narrowband sound centered around the blade passing 
frequency and its harmonics. The modulation spectrum appears in the sidebands of these 
frequencies. The primary subject of this study is tone noise and therefore only the sound 
pressure field within a narrowband is considered. The width (2B) of this band is assumed 
to be small in comparison with the blade passing frequency. It is further assumed that the 
sound pressure field of the (m,n)^^ mode of all the frequencies within such a band can be 
represented by a single mode at the center frequency. The equation for the amplitude of 
such a narrowband mode has been shown in the preceding section, where it also has been 
pointed out that this amplitude is modulated. An eddy convected through a fan creates a 
set of pulsed modes. For a given fan and fan operating condition, these pulses are a function 
of the eddy parameters 

R, 4> eddy center coordinates 
a, £, C| eddy scales 

The results of this section are presented in the form of sound power levels, which are a 
function of the mean squared mode amplitudes. The major task of the analysis is therefore 
to compute the mean squared mode amplitudes resulting from a homogeneous distribution 
of eddies convecting through a fan. 

The mode amplitude as defined by equation (5.1.05) contains only one time dependent term 


1 1 
M 


This function describes the mode amplitude pulse and becomes a Gaussian for large temporal 
scales C/M and large bandwidth B. The determination of the mean square mode amplitudes 
is based on the turbulence representation used in Section 4.3. First, a set of identical eddies, 
Poisson distributed along a streamline is considered. The average occurrence rate at a 
stationary point is Q R d4>dR. Based on Lee (Ref. 36), the mean square value for such a 
pulse distribution is equal to the product of the pulse energy and the average occurrence 
rate. The energy of the mode amplitude pulse is 
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The mean square mode amplitude due to Poisson distributed, identical eddies on a single 
streamline is therefore 
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(5.2.04) 


For a homogeneous field of identical eddies Poisson distributed along all streamlines in 
the annulus and with independent distributions on alt streamlines, the mean square mode 
amplitude becomes 
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(5.2.05) 


In this formulation the interference between the solid annulus walls and the convected 
distortion elements has of course been neglected. 

So far, only eddies of identical scales have been considered. In order to represent actual 
turbulence, the homogeneous field is composed of eddies with different length and velocity 
scales. In the general representation of the turbulent inflow field the scales a, C and c are 
assumed to be random variables. The mean square mode amplitude resulting from such a 
field is 
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with 

P(a, C, c) = joint probability density distribution for the eddy scales 


^mnk 
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= single scale mean square mode amplitude 


In order to make the above integrations feasible, it is assumed that the eddy scales are statis- 
tically independent random variables, each with its own probability density distribution 

P(a) P(£) P(c) 


In order to conform with the model described in Section 4.3, Beta probability density 
functions are used for all three variables 
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(5.2.07) 


with 

t = normalized random variable 
r = gamma function 

and the velocity scale and the length scales are related through the following deterministic 
relationship. 




(5.2.08) 


Based on the assumption of statistical independence for the three random variables the 
integrations can be separated and the equation for the mean square mode amplitude is 
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In Section 4.3 formulas have been shown for the turbulence variance of a homogeneous 
random vortex and eddy field. Based on these equations ((4.3.15) and (4.3.25)), the mean 
square mode amplitudes can be related to the turbulence variance. 
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= turbulence intensity 
= 1. for Gaussian eddy 

= 2/e for Vortex 
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The integrals in equation (5.2.10) are solved by numerical methods. 

The mean acoustic power resulting from a homogeneous random eddy field convected 
througli a fan rotor is computed with equation (5.1.09) based on the mean square mode 
amplitudes. The cutoff criterion is again used to select the propagating modes. It is assumed 
that the cutoff condition at the blade passing frequency and its harmonics is approximately 
valid for all frequencies within the band of width 2B. 

5.2.2 RESULTS 

A parametric study of rotor tone noise due to convected turbulence has been conducted 
based on the model described in the preceding section. The purpose of this study is to 
determine the sensitivity of the sound power level at blade passing frequency to various 
turbulence and fan rotor parameters. The study is based on the fan geometry described in 
Section 5.1.2. The fan operating conditions are again limited to the approach power range, 
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high subsonic tip speeds. The parameters investigated, their ranges and the reference case 
are shown in the following table. 


Parameter 

Range 

Base case 

Number of rotor blades 

12 

46 

23 

Mechanical tip Mach number 

.7 

.874 

.874 

Relative tip Mach number 

.745 

.930 

.930 

Transverse integral length scale (duct radii) 

.01 

.5 

.5 

Streamwise integral length scale (duct radii) 

.1 

100. 

100. 

Distortion velocity component 
Distortion element 

streamwise, 

transverse 

angular 

vortex 

Gaussian eddy 

transverse 

angular 

Gaussian eddy 


The parametric study of rotor tone noise due to convected turbulence has provided the 
following results: 


Turbulence Variance 

Equation (5.2.10) indicates that the sound power is proportional to the relative turbulence 
variance (a/M)2, which is the same for all results shown in this section. 

Transverse Integral Scale 

The effect of the transverse integral scale on the sound power level seems to be rather small, 
as indicated in Figure 37. Especially at large integral scales, the power level decreases at a 
small rate. At first this seems to be in contradiction with the unsteady distortion results 
shown in Section 5.1.2, where the power level dropped off very rapidly with increasing 
scale. It must be emphasized that we are now considering the sound pressure field due to 
a field of distortion elements of various sizes, whereas the results of Section 5. 1 .2 are based 
on single distortion elements. The results of this section simply indicate that even for large 
integral scales there is a sufficient amount of energy in the eddies with smaller scales that 
dominate the contributions to the sound pressure field. Based on the spectral distribution 
of the turbulent energy, the sound power assumes an exponential relationship with the 
integral scale at scales that are sufficiently larger than the scale at the peak power level. The 
exponent is (-2/3), due to the (-5/3) exponent for the inertial subrange of the turbulence 
and the assumption of a constant turbulence variance (see Fig. 10). The peak power level 
occurs at a transverse integral scale similar to the transverse eddy scale predicted by the 
unsteady distortion results. 

Turbulence Components 

In Section 5. 1.2 it has been pointed out that the streamwise turbulence component has a 
larger blade normal velocity component than the transverse turbulence component in the 
rotor tip region. For turbulence of the same variance and spectral distribution, this results 
in higher sound power levels for the streamwise component. For a homogeneous distribution 
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Field of 
Gaussian 
eddies 


Transverse integral scale in rotor blade tip spacings 


Figure 37. Sound Power Level for Homogeneous Distributions of 
Various Distortion Elements 



Transverse integral scale in rotor blade tip spacings 


Figure 38.— Sound Power Level as a Function of Transverse Integral Scale 
for Various Transverse Scale Distributions 
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of Gaussian eddies representing turbulence of the same intensity for the streamwise and the 
transverse turbulence components, the difference in the sound power level is 6.8 dB 
(Fig. 37). 

Distortion Elements, Gaussian Eddy or Vortex 

The sound power level due to a Gaussian eddy field and a vortex field are compared in 
Figure 37. As in the unsteady distortion results, the power levels for the vortex field are 
significantly lower. The reason for these lower levels is the fact that the sound pressure field 
due to a vortex does not couple well into the lower radial order modes. This is demonstrated 
by the modal power distributions shown in Figures 39 and 40. The power levels of the 
highest radial order (n = 3) modes are very similar, but for the lowest radial order (n = 1 ) 
the power levels are much higher for the Gaussian eddy field. This difference decreases with 
increasing power contribution from higher radial order modes, for example, for rotors with 
large numbers of blades and high fan tip Mach numbers. 

Probability Density Distribution (PDD) for the Transverse Scale 

In Section 4.3.3 it has been shown that the PDD for the scales have a significant effect on 
the turbulence spectra. It is therefore of interest to know their effect on the corresponding 
sound power levels which result from the turbulence convecting througli a fan rotor. Three 
FDD’s are considered; the standard PDD (mg = .0, n^ = .0) which has been used to compute 
the bulk of the data presented in this report; a right-sided PDD (m^ = 50., n^ = .0) with its 
highest probability densities close to the maximum scale, and a left-sided PDD (m^ = .0, 
n^ = 50.) with its highest probability densities in the smallest scales. Sound power levels 
for the three distributions as a function of transverse integral scale are shown in Figure 38. 
The largest deviations from the power levels based on the standard PDD occur for the right- 
sided PDD. The latter approximates a delta function at a scale close to the maximum scale, 
and the spectrum for the transverse wavenumber component approaches the shape of a 
Gaussian distribution function. This explains the rapid decrease of the sound power levels 
at large integral scales for this distribution, similar to the results reported in Reference 38. 

Rotor Blade Spacing 

This comparison is based on constant rotor solidity. Sound power levels for fan rotors with 
different numbers of blades are shown in Figure 41 as a function of the transverse integral 
scale. The peak power levels are nearly the same but occur at different transverse integral 
scales. The scales corresponding to the peak power levels are in the order of 25% of the 
rotor blade gap. The following table summarizes the results for the JT9D and JTl 5D engines. 


Engine 

Inlet 

radius 

Ro(cm) 

Number 
of rotor 
blades 

N 

Transverse integral scale at peak power level 
L L N 


Ro 

27tR„ 

Hcm) 

JT9D 

118. 

46 

.035 

.25 

4.1 

JT15D 

27. 

28 

.05 

.25 

1.3 
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Rotor tip mach number = .874 
number of rotor blades 12 


Radial order 
n 

♦ Sum over all radial orders 

• 1 
■ 2 

A 3 



-10 0 10 
Spinning order {m) 


Figure 39.— Modal Power Distribution for a Field of Gaussian Eddies 
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Rotor tip mach number = .874 
number of rotor blades 1 2 

Transverse integral scale = .5 duct radii 


Radial order 



-10 0 10 
Spinning order (m) 


Figure 40.— Modal Power Distribution for a Vortex Field 
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Constant rotor solidity 
Constant turbulence variance 



TRANSVERSE INTEGRAL SCALE IN DUCT RADII 



TRANSVERSE INTEGRAL SCALE IN ROTOR BLADE TIP SPACINGS 


Figure 41. —Sound Power Level as a Function of Transverse Integral Scale 
for Various Numbers of Rotor Blades 


In Figure 42 the cumulative power level is shown as a function of the upper bound on the 
scale distribution. In these results, only a limited segment of the total scale distribution 
has been considered, namely the scales smaller than the transverse scale limit. The cumulative 
power level represents the sound power resulting from ail the eddies smaller than the scale 
limit within a fixed scale distribution. The minus 3 dB point (the half power point) occurs 
at scales in the order of 25% of the rotor blade gap. The results also show that only a 
very small amount of sound power is generated by scales larger than about four times the 
scale at the peak power level. This indicates again that even though the power level is only 
a weak function of the integral transverse scale, the bulk of the sound power is generated 
by a narrow range of transverse scales centered around a transverse scale of about 25% of 
the rotor blade spacing in the fan tip region. 


Tip Mach Number 

Results for three different tip Mach numbers are presented in Figure 43. The results show 
that the peak power level shifts towards larger integral scales with increasing tip Mach 
number. The increase in the power level with the tip Mach number is scale dependent. 

Based on the data the exponent of the Mach number in a power law relationship with the 
sound power would vary between 4 and 7. 

Higher Order Harmonics of the Blade Passing Frequency 

Power levels for the second and fourth harmonic of the blade passing frequency are shown 
in Figure 44. The data demonstrate a clear dominance of the blade passing frequency over 
its harmonics. The presented results are based on a rotor with 12 blades. For the larger 
blade numbers (23 or 46) the blade passing frequency would dominate throughout the whole 
range of investigated transverse integral scales. 


Streamwise Integral Scale 


The effect of the streamwise scale distribution on the mean square mode amplitudes and 
the sound power is represented by the integral over the normalized streamwise scale parameter 
t in equation (5.2. 10). For a given scale distribution, a set of exponents mg and ng, this 
integral is only a function of the parameter combination (B . ®MAX/^)- Therefore the 
mean square mode amplitudes as well as the sound pressure level resulting from a given scale 
distribution can be presented as a function of this parameter combination. The maximum 
scale can be substituted with the integral scale L and the parameter group can be 
presented in various forms. 


with 

T 

S 

^o 

N 



M 




J_ _B ^tip 
27T COq M 


NL 


temporal integral scales 

number of rotor blades cutting through an eddy of length L 
angular blade passing frequency 
number of rotor blades 
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Upper transverse scale limit in duct radii 


Figure 42.— Sound Power Level as a Function of the Upper Limit 
on the Transverse Scale Range 
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Figure 43.— Sound Power Level as a Function of Transverse Integral 
Scale and Fan Tip Mach Number 
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Figure 44.— Sound Power Level as a Function of Transverse Integral 
Scale for Several Blade Passing Frequency Harmonics 
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All parameters are non-dimensionalized with the inlet radius and the speed of sound. For 
the selected scale distribution (mg = ng = .0) the sound power level is presented as a function 
of Q in Figure 45. The nature of the relationship can be explained by the characteristics 
of the error function. For small values of Q the sound power level is proportional to it. 

At large values of Q the sound power level becomes independent of it. The minus 2 dB 
point which represents about two-thirds of the maximum sound power occurs roughly at 
the point where the two limit solutions intersect and the corresponding value for Q is 
about .25. This value can be used to define the conditions for which most of the power of 
rotor tone noise due to convected turbulence is concentrated within a band of width 2 B 
around the blade passing frequency and its harmonics. For the selected fan geometry the 
conditions are: 


7T M 

L> 

2 B 


7T 

T> — 
2B 




S> .25- 


7T ^O M 

L> . 

2 BN Mtip 


If the bandwidth (2 B) is 1% of the blade passing frequency, then the average eddy of integral 
scale L should be cut by 50 rotor blades. For a fixed bandwidth the length scale require- 
ments are independent of the number of rotor blades and the rotor size. 

5.3 ROTOR TONE NOISE DUE TO ATMOSPHERIC TURBULENCE 
5.3.1 BACKGROUND 

It has long been known that steady and unsteady distortions in the inflow field of the fan of 
a stationary engine can cause significant fan noise levels in the vicinity of the blade passing 
frequency and its harmonics (Refs. 1 through 6). These distortions have been attributed to 
the atmospheric turbulence, the ground vortex and the flow field around the test stand 
structure (Refs. 5 through 8). The interference of these distortions with the inlet boundary 
layer and the resulting distortions have also been suggested as a significant cause of fan tone 
noise (Refs. 6 and 20). Most of these distortions do not exist in flight and as a result of this, 
the fan sound pressure field changes significantly between the two conditions. The distortions 
associated with the atmospheric turbulence are present throughout the atmospheric boundary 
layer and therefore could possibly affect fan noise along the approach flight path. Very little 
work has been done in the assessment of fan tone noise due to atmospheric turbulence in 
flight. A Boeing report used simple turbulence and acoustic models to estimate fan tone 
noise due to atmospheric turbulence and showed significant differences between static and 
flight conditions. In this section these differences are evaluated based on the models 
described in the preceding sections. The tone noise generation due to atmospheric turbulence 
convected through a fan is a rather complex process. Many simplifying assumptions had to 
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be made in order to provide a prediction based on analytical models. Due to these simplifica- 
tions, the results must be considered rough estimates at best. The study is only concerned 
with the changes in the sound source, differences in the sound propagation and radiation are 
not considered. 

5.3.2 APPROACH 

In this study it is assumed that the turbulence convecting through a fan evolves from the 
turbulence in the atmospheric boundary layer. In this process a small segment of the entire 
boundary layer is strained by the fan inflow field and is ingested into the fan. The difference 
in the resulting fan tone noise between static and flight conditions is due to two principal 
causes. 


Differences in the inflow contraction between the ambient and fan face conditions. 


• Changes in the atmospheric turbulence with altitude and the surrounding environment. 

In the following outline of the approach, an approximate relationship between the turbulence 
involved in the tone noise generation and the turbulence of the atmospheric boundary layer 
is derived. The results presented at the end of this section are based on the exact equations 
of the models described previously. 

Atmospheric Turbulence 

The description of the atmospheric turbulence is based on the model of Section 4. 1 . The 
turbulence is characterized by its large scale. All integral scales increase with altitude. The 
streamwise integral scale of the vertical turbulence component is 


L^ = .5z (5.3.01) 

The turbulence standard deviations are proportional to the friction velocity u^^, which is a 
function of the wind velocity Uj-gf and the surface roughness Zq. It can be approximated by 


“♦o 



(5.3.02) 


For neutral stability the relationship between the standard deviation of the vertical turbulence 
component and the friction velocity is 

a^=1.3.u*Q (5.3.03) 
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In the first order approximation is independent of altitude. These three relationships 
indicate that the spectral distribution of the turbulence must change with altitude. Due to 
this the turbulence is assumed homogeneous only in horizontal planes. 


In neutral atmospheric conditions, the mean wind velocity can be approximated by 


U = 2.5 u*Q 



“ ^REF 



(5.3.04) 


The preceding formulas are used to describe the atmospheric turbulence in the approximate 
relationship between the turbulence involved in fan tone noise generation and the 
atmospheric turbulence. 

The Inflow Streamtube 

The inflow field of a fan can be described as a streamtube in the ambient flow field. The 
size (radius of this streamtube is a function of the contraction ratio and the radius 
of the inlet Rq 


Ra=RonAT (5.3.04) 

The contraction ratio is defined as the ratio between the mean flow velocity at the fan face 
and the relative velocity between the inlet and the ambient air. The latter is the local wind 
velocity U(z) in static conditions, and the airplane velocity V^p in flight conditions. This 
streamtube is approximately circular in flight and its diameter is in the order of the inlet 
diameter. In static tests, where the wind velocities are required to be small, the streamtube 
becomes very large and thus the mean flow velocity distribution within the tube cannot be 
assumed uniform. 

For several reasons, it is desireable to define the turbulence associated with distortion elements 
of transverse scales smaller than the radius of the inflow streamtube. 

• In static tests, elements with scales in the order of the streamtube radius are only 
partially affected by the inflow contraction. The increase in the energy of the 
transverse turbulence component of such elements will not be as large as for smaller 
elements, completely contained within the streamtube. 

• In the turbulence model of Section 4.3 which forms the basis of the rotor tone noise 
model of Section 5.2, it is assumed that the streamwise and transverse length scales 
are statistically independent random variables. This assumption can lead to unrealistic 
results, if the range of scales contributing significantly to rotor tone noise comprises 
only a small segment of the total distribution of scales. 
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• Distortion elements of the transverse turbulence component with scales in the order 
of their separation distance from the inlet surface are suppressed. 

• The results of Sections 5.1 and 5.2 indicate that the distortion elements with transverse 
length scales in the order of one-half of the duct radius or larger contribute only a 
negligible amount to the sound pressure field resulting from convected turbulence. 

For these reasons, the study of rotor tone noise due to atmospheric turbulence will be 
based on the turbulence resulting from the scales smaller than one-half of the streamtube 
radius. 


a < .5 

The atmospheric turbulence due to a limited scale range is determined with the inflow 
turbulence model of Section 4.3. Spectra based on this model representing different 
ranges of scales are shown in Figure 3 1 . Based on this simple model, the variances 
and the integral scales L of the limited (subscript A) and the full range of scales are related 
by the following equation: 



This relationship is of course based on spectra with an inertial subrange exponent of 
(-5/3). 

Various investigators (Refs. 26 and 27) have shown that atmospheric turbulence becomes 
isotropic at wavelengths small in comparison with the altitude above ground. Teuniessen 
(Ref. 2 1 ) suggests that isotropy occurs for wavelengths < z/3. 

In the flight condition, the mean flow velocity is uniform throughout the streamtube and 
the maximum scale of the limited scale range is only a small fraction of the altitude. The 
turbulence based on the limited scale range is therefore homogeneous and isotropic. 

In static conditions, the situation is more complex. The wind velocity varies significantly 
across the streamtube. Depending on the streamtube size and the engine height above the 
ground, the maximum scale of the limited scale range can be in the order of the altitude 
above the ground. The precontraction turbulence is therefore neither homogeneous nor 
isotropic. The flow contraction model is, however, based on isotropic, homogeneous 
turbulence. Some simplifying assumptions must be made in order to conform with the 
requirements of this model. 

At low altitudes, the vertical turbulence component has a smaller variance and integral scale 
than the horizontal components. In the spectra, this difference manifests itself in the 
dominance of the horizontal components at small wavenumbers. It is assumed that this 
large scale horizontal motion does not contribute to the turbulence in the fan inflow field. 
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At large wavenumbers the spectra of the three components merge. A limited segment of the 
atmospheric turbulence is assumed to be isotropic. The integral scales and variances of this 
turbulence are assumed to be equal to the corresponding values for the vertical component. 

L=L^ a = a^ (5.3.06) 

The range of scales that is used to represent the turbulence in the fan inflow field is a 
function of the relative size of the streamtube and the integral scale L^. For the selected 
scale distributions with the exponents m and n set equal to zero, the ratio between the 
maximum scale and the integral scale is 1.41. Therefore, if the maximum scale of the considered 
scale range is one-half of the streamtube radius, 

^max “ ^A 

Then the integral scale of the precontraction turbulence is 

— .355 • 

If this integral length scale is smaller than the integral length scale of the vertical turbulence 
component then only a limited scale range of the actual turbulence is considered. But if 
the length scale based on the above equation is larger than L^, then the full spectrum of 
the vertical turbulence component is used. These two conditions are described by the 
following equations. 


L^> .355 



'.355 Ri\ 

— 1 ,La = .355 Ra 

. ^W / 


L^<.355 


a A = CT, 


, (5.3.07) 


Based on the approximate equations for the atmospheric turbulence and the assumption 
> .355 R^, the standard deviation for the fan inflow turbulence is 


^A ^w 



1/3 



(5.3.08) 
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(5.3.09) 
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Inflow Contraction 


The effects of the inflow contraction are estimated based on the contraction models of 
Section 4.2. The model is used to compute the post contraction integral scales and 
variances. In the first order approximation, it is assumed that the scales are changed propor- 
tional to the corresponding contraction ratios. The post-contraction standard deviations of 
the turbulence are approximated with PrandtTs formulas 

^2B"v^-ota (5.3.10) 

which lead to the following approximate formulas for the fan face turbulence intensities. 



'2B 


Ui 




(5.3.12) 


The above equations represent first order approximations for the static test conditions. 
In-flight conditions the airplane velocity enters the contraction ratio 


U. 
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(5.3.13) 


and the corresponding approximate equations for the turbulence standard deviations are: 
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(5.3.14) 
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(5.3.15) 
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Fan Tone Noise 

Differences in fan tone noise due to atmospheric turbulence between static and flight condi- 
tions are evaluated with the model described in Section 5.2. In this comparison the turbulence 
intensities, the streamwise and transverse length scales and the general character of the spectral 
distributions are taken into account. 

In static conditions the turbulence at the fan face is dominated by the transverse component. 
Due to the large contraction in the inflow field the turbulence is well represented by a vortex 
field. In flight conditions all turbulence components are of similar magnitude due to the 
small inflow contraction ratio. According to the results of Section 5.2 the streamwise 
turbulence component causes significantly higher noise levels in such conditions, due to the 
fact that it has a larger blade normal velocity component. For the small contraction ratios 
typical for the flight condition the turbulence is represented with distributions of Gaussian 
eddies. 

5.3.3 RESULTS 

A limited parametric study of rotor tone noise due to atmospheric turbulence has been 
conducted to estimate the difference between static and flight conditions. The study is 
based on the exact equations of the models described in the preceding sections and on the 
approach outlined in Section 5.3.2. It is limited to the approach power range and the rotor 
geometry is described in Section 5.1.2. In Section 5.2.2 rotor tone noise due to homogeneous 
but anisotropic turbulence has been investigated for different fan rotors and operating 
conditions. The following study is primarily concerned with the effects of the parameters 
that describe the atmospheric turbulence. It is based on the fan size of a JT9D engine and 
is centered around two reference conditions representing static and flight tests. 


Reference conditions 

Static 

Flight 

Altitude 

4.88 m 

120 m 

Reference altitude 

4.88 m 

10m 

Wind velocity at reference altitude 

2.68 m/s 

5.14 m/s 

Richardson number at reference altitude 

.0 

.0 

Roughness scale 

.2 m 

.02 m 

Boundary layer thickness 

600 m 

Equivalence altitude 

600 m 

Airplane velocity 

.0 

82.3 m/s 

Rotor radius 

1.16m 

Axial Mach number at fan face 

.35 

Relative tip Mach number 

.93 

Number of rotor blades 

23 


Bandwidth 2B 

1%of BPF 
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The static condition represents typical test stand operation. The engine axis is 4.88 m 
(16 ft) above the ground. The wind velocity at engine axis heiglit is 2.68 m/s (6 MPH) 
which represents an upper limit on winds from the forward arc. This limit is usually lower 
for wind directions at an angle relative to the inlet axis and drops further for crosswinds and 
winds from the aft arc. The flight condition is based on landing approach at the FAR Part 36 
measuring station. The wind velocity corresponds to the upper limit of the test window. 

According to Reference 22 a roughness scale of .02 meters is typical for airport environments 
and is used for the flight condition. A value of .2 meters has been selected for the static 
test stand environment. 

A boundary layer thickness of 600 meters as suggested in References 21 and 23 has been 
used. Throughout this study it is assumed to be invariant and the same for both conditions. 

Since this study is concerned with narrowband noise only, the rotor noise within a small 
band centered around BPF is considered. A bandwidth of 1% of BPF has been assumed. 

Comparison of the Two Reference Conditions 

The difference in fan narrowband noise due to atmopheric turbulence between static and 
flight conditions is due to the following discrepancies: 

• At the low wind velocities of static tests (2.68 m/s, 6 MPH) a large flow field contrac- 
tion exists between the ambient and the fan face conditions. This contraction causes 
an increase in the transverse turbulence velocities. The corresponding flow contraction 
in flight is small due to the airplane forward velocity. 

• Due to the difference in the contraction ratios, the turbulence originates from larger 
scales and therefore higher turbulence energy levels in static tests. 

• Due to the difference in the contraction ratios, the streamwise integral scale is much 
larger in static tests. Therefore a larger amount of the turbulence energy contributes 

to narrowband sound in a small band centered around the BPF. In flight the streamwise 
integral scale is small and therefore only a fraction of the total turbulent energy con- 
tributes to narrowband sound at BPF. 

• Due to the difference in contraction ratios, the streamwise turbulence component 
dominates the noise generation in flight, whereas the transverse component dominates 
in static conditions. 

• Within the surface layer of a neutral atmospheric boundary layer the turbulence 
variances are nearly independent of height, but the integral length scales increase with 
altitude. At large wavenumbers this results in a decrease of the energy within a given 
wavenumber band with altitude. In other words, the energy associated with a turbu- 
lence scale decreases with altitude. 
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All these differences except the fourth one lead towards higher fan tone noise levels in 
static tests. 

The approximate intensities of the turbulence involved in the tone noise generation can be 
estimated with equations (5.3.12) and (5.3.14). Based on the assumptions that Z = Zj^gp 
in static conditions the ratio between the turbulence intensities for flight and static 
conditions is 



with subscripts F = flight 
S = static 

For the two reference conditions this ratio is 



(5.3.16) 


The difference in the narrowband sound power level W at BPF between static and flight 
conditions is usually dominated by the change in the intensity of the turbulence involved 
in the tone noise generation. The power difference between the two reference conditions 
is 


W = 27.3 dB 

This difference is representative for the changes in narrowband fan noise due to atmospheric 
turbulence between static and flight conditions. 

The effect of the various independent variables on the sound power level are discussed in 
the following paragraphs. The sound power level is normalized with its value for the static 
reference condition. 

Altitude 

In flight, the intensity of the turbulence involved in tone noise generation decreases due to 
the increase in the integral scales and the slow decrease in the turbulent shear stresses 
(Figure 46). The corresponding change in the sound power level is shown in Figure 47. 
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Figure 46.— Estimated Turbulence Intensities at the Fan Face Based on 

Scales Involved in Rotor Tone Noise Generation as a Function 
of Altitude Above the Ground 
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Figure 47.— Predicted Narrowband Sound Power Level W at BPF Due to 

Atmospheric Turbulence as a Function of Altitude Above Ground 
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In static conditions the decrease in the power level with altitude is even stronger. This is 
due to the decrease in the contraction ratio with increasing altitude. This is an indication 
that static test conditions can be improved by increased engine height above the ground 
level. 

In static conditions the fan inflow field is not uniform due to the variation of the mean wind 
velocity with altitude, which results in a variation of the contraction ratio throughout the 
inflow field. This variation in the contraction ratio coupled with the variation in the 
turbulence properties with altitude result in a non-homogeneous turbulence field at the fan 
face. The models used in the present study are not sophisticated enough to evaluate such 
complex flow fields. But by assuming local uniformity in the mean flow field and homo- 
genity in the turbulence field, the turbulence intensity at the fan face can be estimated for 
streamlines originating at different levels above the ground (Figure 46). The results indicate 
a significant increase in the intensity of the transverse turbulence component for flow 
segments originating closer to the ground. The change in the slope of the curves occurs at 
the point where the integral length scale of the considered scale range is equal to the integral 
scale of the vertical turbulence component. The corresponding sound power levels 
(Fig. 47) indicate increasing levels with decreasing engine axis height above the ground. In 
an actual static test situation with an inhomogeneous inflow field, the highest contributions 
to the sound pressure field are therefore expected to come from the bottom segment of the 
inlet. Rotor blade mounted transducer data reported by Hanson (Ref. 39) do indeed show 
a variation of the standard deviation of the pressure signals along the circumference of the 
fan face, with the highest values in the bottom segment of the inlet. 

In most static test conditions, the inflow streamtube touches the ground. In that case, a 
stagnation point forms on the ground and along the streamline from this point the contrac- 
tion ratio becomes infinite. The vorticity associated with the horizontal turbulence compo- 
nents, which have an intensity twice as large as the vertical component close to the ground, 
are very effectively amplified by the transverse contraction around such a stagnation point 
streamline. The contraction of the atmospheric turbulence in the flow field around such 
streamlines leads to a phenomena frequently referred to as a ground vortex. It is very likely 
that the highest turbulence intensities at the fan face occur in the vicinity of such stagnation 
point streamlines. The models used in this study do however not allow the prediction of the 
resulting intensities. Similar stagnation point situations can occur on surfaces of the test 
stand structure in the vicinity of the inlet. 

Surface Roughness Scale 

The surface roughness scale affects the mean velocity profile as well as the standard deviations 
of the turbulence velocity components. Its effect on the length scales of the ingested turbu- 
lence can be neglected and its effect on the intensity of the ingested turbulence as well as 
on the sound power can be estimated with the approximate equations (5.3.12) and (5.3.14). 


for Z = Z 
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The ratio Zj^£p/Zq is usually quite large and therefore the effect of Zq is relatively small. 

In the flight condition, Zq is reasonably well defined. Most landing approaches occur in 
headwinds and therefore the boundary layer evolves along the airfield. Suggested values for 
the roughness scale of airports range from .02 to .05 meters. 

In static conditions it is more difficult to define a surface roughness scale. In most static 
test stand environments the surface cover changes abruptly in the vicinity of the test arena. 
Such a change will initiate a new internal boundary layer, a situation that exceeds the 
limitations of the simple atmospheric boundary layer model. Roughness elements with scales 
in the order of the engine height above the ground (buildings, test stands, trees, forests) 
are present in the vicinity of most test stands. The flow around such objects creates large 
scale, high intensity turbulence. The validity of the simple atmospheric boundary layer 
model is limited to altitudes higher than the size of the roughness elements. Usually the 
surface roughness varies around the test stand and therefore the mean velocity as well as the 
turbulence intensity change with the wind direction. These problems make it difficult 
to predict atmospheric turbulence in static test stand environments. 

Inflow turbulence was measured at the fan face of a JT9D-7 engine during two static tests 
which form part of a joint P&WA and Boeing noise reduction program. The two tests were 
conducted at different test facilities. The turbulence was measured with split film and 
X-wire anemometers located at several radial positions on a rake in the bottom segment of 
the inlet. The statistical averages of the turbulence measured at the fan face, averaged over 
several conditions and all radial locations are listed in the following table. 
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The turbulence intensities are presented in Figure 46 together with results predicted by the 
analytical model for various surface roughness scales. The results indicate that the predicted 
and the measured turbulence intensities for the transverse component are of the same order 
of magnitude. An extremely large discrepancy exists however between the predicted and 
measured intensities of the streamwise turbulence component. This difference can be 
partially explained with the high noise floor in the streamwise turbulence component signal 
measured by the split film anemometers, which affected the higher of the two levels. It 
also could be due to density fluctuations in the highly anisotropic turbulence. However, 
many investigators have reported, that the attenuation of the streamwise turbulence component 
is not nearly as large as predicted by the analytical flow contraction models. Nonetheless, 
the separation between the two turbulence component intensities is large enough to make 
the transverse component the dominant one in the fan tone noise generation in static tests. 


104 




The narrowband sound power at BPF due to atmospheric turbulence in static tests has been 
computed based on the measured turbulence length scales and the intensities of the transverse 
turbulence component. The power is computed with the assumption of a homogeneous 
turbulence distribution. Since it is based on the turbulence characteristics measured in the 
bottom of the inlet, it overestimates the sound power level to be expected in static test 
conditions. The results are presented in Figure 47 together with the sound power levels 
predicted by the analytical model. The two sets of data agree reasonably well. Based on 
these results, it is concluded that the surface roughness scale in the test stand environment 
should be minimized. Large roughness elements (buildings) should be moved away from the 
test arena. 

Mean Wind Velocity 

The reference wind velocity affects the mean wind velocity as well as the standard deviations 
of the turbulence components. Its effect on the integral scales of the turbulence involved 
in the generation of narrowband fan noise is minimal. Its effect on the corresponding 
turbulence intensities can be estimated with equations (5.3.12) and (5.3.14). 



The results based on the exact equations are presented in Figure 48. The corresponding 
sound power levels are presented in Figure 49. The results indicate that in-flight power 
levels increase rapidly with the mean wind velocity, but in all situations they remain substan- 
tially lower than the ones occurring in static tests. In static test conditions the effect of the mean 
wind velocity is less significant. But lower wind velocities provide lower turbulence intensities 
at the fan face and are therefore preferred in static test conditions. 

Atmospheric Stability 

In an atmosphere with an adiabatic lapse rate, in neutral stability, the turbulence is sustained 
by the shear in the mean flow. The associated energy transfer from the mean flow to the 
turbulence is proportional to the velocity gradient and therefore decreases with altitude. 

In stable or unstable conditions energy is subtracted or added to the turbulence through 
buoyancy forces. This energy is assumed to be invariant with altitude. The absolute value 
of the flux Richardson number, the ratio between the turbulence production by buoyancy 
and by shear increases therefore with altitude. 

This is reflected in the predicted intensities for the turbulence involved in fan tone noise 
generation shown in Figure 50. They are based on a mean wind velocity of 4.12 m/s 
(8 kts) at 6.1 meters (20 ft) altitude and the range of Richardson numbers has been selected 
based on data presented in Reference 16, which indicate that at the selected altitude and 
mean wind velocity Richardson numbers smaller than -.5 have an occurrence probability 
of less than 10%. The same set of data indicates that the average Richardson number is 
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Figure 48.— Estimated Turbulence Intensities at the Fan Face Based on 

Scales Involved in Rotor Tone Noise Generation as a Function 
of the Wind Velocity at 10 Meters Above the Ground 
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Figure 49.— Predicted Narrowband Sound Power Level W at BPF Due to Atmospheric 
Turbulence as a Function of the Wind Velocity at 10 Meters 
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slightly positive. The results shown in Figure 50 demonstrate the negligible effect of thermal 
stability on the static conditions and its substantial effect on the flight conditions. The 
corresponding sound power levels are presented in Figure 5 1 . The gap between static and 
flight conditions is substantially reduced in unstable conditions, but is still in the order of 
20 to 25 dB. 

Airplane Mach Number 

The analytical model has been used to estimate the effects of airplane forward velocity on 
the ingested turbulence as well as on the associated fan narrowband sound power at BPF. 

With increasing forward velocity, the inflow contraction ratio is reduced and therefore the 
contraction effects on the turbulence decrease. At the same time the size of the inflow 
streamtube is reduced which in turn reduces the scale and therefore the energy of the pre- 
contraction turbulence. The resulting intensities for the turbulence involved in the fan 
narrowband noise generation are presented in Figure 52. They are based on the size and 
location of a JT9D engine installed in a 747 aircraft. The reference wind velocity was 
assumed 2.68 m/s (6 MPH) at 4.88 meters (16 ft). The results are based on the conditions 
along a streamline that enters the fan horizontally at engine axis height. They demonstrate 
the significant reduction in the intensity of the transverse turbulence component during the 
initial phase of the ground roll (M^p < .05). They also show insignificant changes for the 
range of aircraft velocities during landing approach. The corresponding narrowband sound 
power levels at BPF are shown in Figure 53. They demonstrate the significant sound power 
reduction (20 dB) during the initial aircraft acceleration to a Mach number of .05 and the 
negligible changes over the range of airplane approach velocities. Significant improvements in 
the fan inflow field can therefore be achieved even at small vehicle velocities (M^.^ > 05). 

^ dp 

Fan Size 

In the evaluation of fan tone noise due to atmospheric turbulence of different scale fans, 
the results cannot simply be scaled. This is due to the fact that several significant scales 
affect the results, the fan radius, the surface roughness scale, the test stand height, the flight 
path altitude, and the thickness of the atmospheric boundary layer. According to the 
models used in the prediction of the narrowband sound due to atmospheric turbulence, the 
acoustic results are primarily a function of the integral scales and the intensities of the 
turbulence at the fan face. For most conditions, a given engine height and contraction ratio, 
the model predicts that the integral length scales are proportional to the fan rotor radius. 
Therefore, the non-dimensionalized integral scales are independent of fan rotor scale. Based 
on this result the sound power could be scaled with the square of the rotor radius. However, 
according to equations (5.3.12) and (5.3.14) the turbulence intensities are a function of the 
fan rotor radius Rq. 
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Figure 51.— Predicted Narrowband Sound Power Level W at BPF Due to 
Atmospheric Turbulence as a Function of Atmospheric 
Stability 
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Figure 53. —Predicted Narrowband Sound Power Level W at BPF Due to 
Atmospheric Turbulence as a Function of Airplane Mach 
Number M^p 
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The predicted turbulence intensities shown in Figure 54 reflect this simple relationship. 

The sound power due to atmospheric turbulence convected through a rotor is proportional 
to the (8/3) power of the fan rotor radius (Fig. 55). In a first-order approximation the 
difference in rotor tone noise between static and flight conditions seems to be independent 
of fan rotor size. Since the sound power due to atmospheric turbulence is proportional to 
the (8/3) power of the rotor scale, it becomes more dominant with respect to other fan tone 
noise sources with increasing fan scale. This change in the separation relative to other sources 
is in the order of 4 dB between the JT9D (Rq = 1.18 m) and the JT15D (R^ = .27 m). 

Rotor Blade Spacing 

Results presented in Section 5.2.2 indicate that changes in the fan tone sound power level 
with the number of rotor blades, or in other words with the rotor blade spacing, are minimal 
for transverse integral scales larger than 3% of the rotor radius. Measured as well as predicted 
transverse integral scales are considerably larger. The predicted results shown in Figure 55 
demonstrate the relatively small effect of the rotor blade number. The level differences 
apparent in the results for the approach flight condition are primarily due to the change in 
the bandwidth which has been assumed proportional to the blade passing frequency. 

Figure 55 also includes predictions for the JT9D and JTl 5D engines. For both engines the 
static to flight differences in rotor tone noise due to atmospheric turbulence are predicted 
to be in the order of 30 dB. 

Fan Tip Mach Number 

In Section 5.2.2 it has been shown that for constant turbulence intensity and scales, the 
fan tone noise power due to convected turbulence increases proportionally with the fourth 
to seventh power of the fan tip Mach number. Based on equations (5.3.12) and (5.3.14) 
the effect of the axial Mach number at the fan face on the turbulence intensities is 



In both conditions the turbulence intensity is reduced with increasing fan face velocity and 
therefore with the fan tip Mach number. This reduces the exponent in the relationship 
between fan tone sound power and fan tip Mach number by 2/3 in the static condition and 
1 1 /3 in the flight condition. Within the approach power range the difference in rotor tone 
noise due to atmospheric turbulence between static and flight conditions is expected to 
increase with fan tip Mach number. 
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Figure 54.— Estimated Turbulence Intensity at the Fan Face Based on 
Scales Involved in Tone Noise Generation as a Function 
of the Fan Radius 




Figure 55 -Narrow Band Sound Power at BPF Due to Atmospheric 
Turbulence as a Function of Fan Rotor Radius R^ and 
Number of Rotor Blades N 



6.0 INLET ANGLE OF ATTACK EFFECTS ON FAN NOISE 


6.1 BACKGROUND 

The flow field in an inlet operated at an angle of attack is asymmetric. The mean flow 
velocity, as well as the boundary layer, vary circumferentially throughout the inlet. The 
resulting circumferential distortions in the flow field at the fan face may affect the sound 
pressure field emitted by the fan. The more significant distortions in the vicinity of the 
inlet surface between the highlight and the throat plane may affect the sound propagation 
from the fan to the far field. In connection with the improvement of test methods for 
static fan noise tests, it is necessary to know the significance of the acoustic effects resulting 
from inlet angle of attack induced flow distortions. If these effects are significant, an 
attempt to simulate them in static conditions might be necessary. 

References 12 and 13 report fan noise results from wind tunnel tests concerned with angle 
of attack effects. The results indicate that inlet operation at an angle of attack will affect 
the directivity and possibly the sound power generated by a fan. Results from an in-flight 
investigation of fan noise reported in Reference 3 indicate that the inlet angle of attack 
might have an effect on fan tone noise at the blade passing frequency and its harmonics. 

During a joint flight test conducted by P&WA and Boeing, a limited amount of information 
concerned with the inlet angle of attack effects on fan noise has been acquired. It is the 
purpose of this section to describe the results of this investigation. 

6.2 TEST DESCRIPTION AND INSTRUMENTATION 

The inlet angle of attack effects on fan noise were investigated during a flight test, which formed 
part of a Joint P&WA and Boeing Noise Reduction Program. The test was conducted in 
Seattle with a 747 aircraft powered by JT9D-7 engines with hardwall nacelles. One of the 
inboard engines was extensively instrumented with static pressure sensors, rotor blade 
mounted transducers and inlet and fan duct microphones. In addition to the standard test 
instrumentation, the airplane was equipped with fuselage-mounted microphones and split 
film anemometers in various locations. The investigation of the inlet angle of attack effects 
is based on aircraft and engine mounted instrumentation only and the evaluation is limited 
to the sound pressure field of the instrumented, inboard engine. In order to minimize the 
effects of parameters other than the inlet angle of attack on the data, the corrected fan speed 
of the inboard engines was kept constant and changes in the aircraft velocity were kept to a 
minimum. The inlet angle of attack changes were obtained by variation of flap settings, 
outboard engine power settings and landing gear position. The investigation is limited to 
two typical power settings for approach and takeoff. The conditions are summarized in the 
following table. 


116 



Condition 

Flaps 

deg 

Landing 

gear 

Airspeed 

m/s 

Inboard 

engine 

NIV^ 

RPM 

Inlet angle 
of attack 

deg 

1 

5 

down 

90. 

3310. 

17.5 

2 

10 

up 

90. 

3300. 

15.5 

3 

20 

up 

91. 

3330. 

11.5 

4 

25 

up 

104. 

3350. 

5.5 

5 

10 

up 

90. 

2425. 

15.5 

6 

20 

up 

88. 

2420. 

12.5 

7 

25 

up 

89. 

2415. 

10.5 

8 

25 

up 

104. 

2425. 

5.5 


During normal takeoff and landing maneuvers the inlet angle of attack varies between 7 and 
15 degrees. The range of angles evaluated during this test therefore covers the normal operating 
conditions. 

The following instrumentation has been used in the evaluation of the angle of attack effects. 
Static Pressure Sensors 

Static pressure ports were located in two planes, upstream and downstream of the fan rotor 
(Fig. 56). The upstream location was separated 9 cm from the fan leading edge and the 
downstream position was 1 5 cm from the trailing edge. Each plane contained 12 equispaced 
pressure ports. Each pressure was measured by an individual pressure sensor and the signals 
were sampled at a rate of 2.5 times per second. The purpose of these measurements was to 
determine the mean flow distortions at the fan face and to observe changes in the static 
pressure distribution downstream of the fan rotor. 

Rotor Blade Mounted Pressure Transducers (BMT) 

Thin film dynamic pressure transducers were mounted on several rotor blades. The diaphragm 
of these sensors is about 1 mm in diameter and the frequency response extends to above 
20 KHz. The transducers were bonded to the pressure side of the blades at various radial 
locations close to the blade leading edge (Fig. 58). The transducer lead wires were run along 
the blade pressure surface to the blade hub and into the nose cone, where they were connected 
to the rotating telemetry system. The signal emitted by the telemetry transmitter was picked 
up by an antenna embedded in the surface of the inlet. From there the signal was carried 
by fixed wiring to a bank of FM receivers, where the signals were separated and passed to the 
proper recording channels. The purpose of these dynamic pressure measurements was to 
detect changes in the steady and unsteady velocity distortions at the fan face. 

Inlet and Fan Duct Dynamic Pressure Transducers 

Dynamic pressure transducers were mounted at various locations in the inlet as well as in 
the fan duct (Figs. 61 and 64). These sensors were used to evaluate changes in the fan 
sound pressure field due to the variation of the inlet angle of attack. 


117 






Fuselage Mounted Microphones 

Twelve .5 inch condenser microphones were mounted in window panes of the test aircraft 
(Fig. 67). They were used to monitor changes in the far field sound pressure distribution 
that could be associated with inlet angle of attack changes. 

All dynamic signals were recorded on tape at 76.2 cm/s. The results presented in the next 
paragraph are based on post test narrowband analysis of these records and represent averages 
for signals of 35 seconds duration. 


6.3 RESULTS 

The results from this test can, at best, provide a partial assessment of the inlet angle of attack 
effects on fan noise. The mean flow field in an inlet operated at an angle of attack is cir- 
cumferentially distorted. The highest velocity distortions occur in the area between the 
highlight and the throat close to the inlet surface. As a result of these distortions in the 
mean flow field corresponding distortions in the propagation must be expected, and this may 
lead to a distortion in the far field acoustic pressures. Very likely such propagation distortions 
cannot be observed with the fuselage mounted microphones. The results of this test will 
therefore at best, provide an assessment of the inlet angle of attack effects on the fan noise 
at its source. 

Static Pressure Sensors 

The circumferential variation of the static pressure at the inlet surface upstream of the fan 
is shown in Figure 57. The static pressure distortion increases with the angle of attack at 
both power settings. The results show a decrease in the static pressures in the bottom of the 
inlet. Results from an earlier flight test show that outside of the thin boundary layer the 
total pressure variations along the inlet circumference are considerably smaller than the 
corresponding static pressure variations. Based on this information, it is concluded that the 
low static pressures in the bottom of the inlet represent a high velocity region. This 
positive streamwise velocity distortion increases with the inlet angle of attack. 

The static pressure distribution downstream of the fan is difficult to interpret without the 
knowledge of the corresponding total pressure distribution. It therefore, is omitted in the 
discussion of the inlet angle of attack effects. 

Blade Mounted Transducers 

The results are based on four transducers located on the same blade at different radial 
locations (Fig. 58). A typical power spectrum of a blade mounted transducer signal is 
shown in Figure 59. The narrow spikes in this spectrum represent the periodic pressure 
fluctuations on the fan blades that result from steady or quasi-steady distortions in the fan 
inflow field. Figure 60 presents the five lowest harmonics of the BMT signals as a 
function of the inlet angle of attack for the four transducers. The data indicate an angle of 
attack related increase in the power of the two lowest harmonics. This increase is greatest 
at the edge of the boundary layer and becomes smaller towards the hub. Changes in the 
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Figure 57.— Circumferential Variation of the Static Pressure at the 
Inlet Surface Upstream of the Fan 
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higher order harmonics do not seem to be related to angle of attack effects. The distortion 
represented by the changes in the two lowest harmonics has a circumferential scale in the 
order of one-half of the circumference and a radial scale in the order of 1 0% of the blade 
span. Such a distortion is either due to the inlet boundary layer or the rotor tip flow. Based 
on the BMT and static pressure data combined, it is concluded that a positive streamwise 
velocity distortion exists close to the inlet surface in the bottom of the inlet. It covers the 
inlet boundary layer as well as part of the mean flow and increases in intensity with the 
inlet angle of attack. At takeoff power an opposite trend is apparent in the top of the inlet. 
These phenomena cannot be explained as inlet boundary layer effects and it is suggested 
that they are due to changes in the rotor blade tip flow. They could result from a change 
in the circumferential distribution of the rotor tip clearance brought about by the asymmetric 
aerodynamic load distribution on the inlet nacelle and the oval shape of the rotor rub strip. 
Relative motion between the fan axis and the fan case under asymmetric aerodynamic loads 
has been observed previously (Ref. 14). 

Inlet and Fan Duct Microphones 

The inlet and the fan duct sound pressure fields were evaluated with flush mounted surface 
transducers only. The recorded data therefore represent only the surface pressure fluctuations 
and their significance in the evaluation of the total sound pressure field should not be over- 
estimated. 

Changes in the broadband noise levels in the inlet as well as in the fan duct were in general 
less than ±1. dB and could not be related to angle of attack changes. This was the case 
for both power settings and the full range of angles of attack. The sound pressure levels at 
the blade passing frequency and its harmonics exhibited large variations between the 
different angle of attack conditions. A typical set of data for the second harmonic at approach 
power in the inlet is presented in Figure 62. This set of data demonstrates the seemingly 
random scatter that does not show any correlation with the inlet angle of attack changes. 

It also shows the large spatial variations in the inlet sound pressure field measured by the 
six coplanar sensors located close to the inlet throat. Such variations result from the 
interference of modes of different spinning order. The scatter in the data could be due to 
changes in the interference pattern between the sound pressure field resulting from the angle 
of attack related distortions and the sound pressure field resulting from a different source. 

It also could be due to the unintentional variation of a parameter other than the angle of 
attack that significantly affects fan tone noise. Experimental data from static tests with 
inflow control structures have demonstrated large, rapid variations in the narrowband sound 
pressure levels with fan rotational speed. It is assumed that small differences as well as 
fluctuations in the fan rotational speed have contributed to the scatter in the data. 

Average narrowband sound pressure levels based on the seven transducers in the fan duct 
and on the six inlet transducers close to the throat are presented in Figures 63 and 65. 

The spatial averaging has eliminated a large portion of the total variation in the signals. 

The variations over the full range of angles of attack is less than 2 dB for most signals, and 
no correlation between the average sound pressure levels and inlet angle of attack is 
apparent. 


122 



K1 1 (98% span) 


K12 (96% span) 




Harmonic order 




I 1 L_ 

5 10 15 


r K13 (91% span) 


K14 (87% span) 




Harmonic order 


1 


2 


4 



3 


i — 1 i_ 

5 10 15 


Inlet angle of attack (deg) 


Figure 60.— Blade Mounted Transducers, Spectrum Levels of the Five Lowest 
Distortion Orders for Sensors at Various Radial Locations 
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Figure 62. —Sound Pressure Level for Second Harmonic of Blade Passing 
Frequency in the Inlet at Approach Power 
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Figure 63.— Average Sound Pressure Levels Based on Six Locations Near the Inlet 
Throat for Harmonics of the Blade Passing Frequency 



An interesting phenomenon was observed at takeoff power in the inlet sound pressure field. 
At small angles of attack, the buzzsaw spectrum is nearly the same at all circumferential 
locations. But with increasing angle of attack, the differences between the spectra in the 
bottom and in the top of the inlet become larger. This is demonstrated by the spectra shown 
in Figure 66. They represent the sound pressure signals of the sensors in the top and in the 
bottom of the inlet at the smallest and the largest angle of attack condition. At the inter- 
mediate locations, on the side of the inlet as well as for the intermediate angle of attack 
conditions, the spectra fall between the extremes. Similar changes in the buzzsaw spectra 
are also observed at the intermediate axial location, between the throat and the fan face. 
Based on this observation, it has been concluded that the asymmetry in the buzzsaw noise 
develops in the vicinity of the fan and is not due to velocity distortions in the throat area. 

It is very likely that this asymmetry is the result of the mean flow distortions in the vicinity 
of the fan face observed with the static pressure sensors and the blade-mounted transducers. 

Fuselage Mounted Microphones 

The signals of the fuselage mounted microphones represent the pressure fluctuations due to 
several different sources; inboard and outboard engine noise, fuselage boundary layer noise, 
landing gear and other aerodynamic noise. The broadband noise floor below 5 KHz is 
dominated by fuselage boundary layer, landing gear and other aerodynamic noise. For this 
reason the evaluation of the inboard engine noise field had to be limited to narrowband 
tones. A typical set of data representing the second harmonic of the blade passing frequency 
at approach power is shown in Figure 68. These data reflect the same scatter problem that 
was observed in the inlet and fan duct pressure signals. No attempt has been made to find 
correlations between the changes in the sound pressure level and the inlet angle of attack 
variations. 

The main results based on this limited investigation are: 

• The most significant changes in the flow field at the fan face associated with changes 
in the inlet angle of attack are very likely due to a change in the rotor blade tip flow, 
rather than a change in the inlet boundary layer. 

• An asymmetry in the buzzsaw sound pressure field that increases with the inlet angle 
of attack has been observed. 

• The effects of inlet angle of attack on broadband noise as well as narrowband sound 

at the blade passage frequency and its second harmonic seem to be in the order of 2 dB 
or less for the sound pressure field at the inlet and fan duct surfaces. 

All conclusions apply only to angles of attack conditions representative of CTOL, subsonic 
aircraft operating conditions. 
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7.0 conclusions: recommendations and summary of results 


Based on a literature search, an atmospheric turbulence model has been selected that 
predicts turbulence intensities, integral scales and spectra as a function of the altitude above 
the ground, the mean wind velocity, the surface roughness and atmospheric stability. An 
existing flow contraction model has been extended to provide all the one-dimensional spectra 
of the post-contraction turbulence for simple, isotropic pre-contraction turbulence. Narrow- 
band fan noise resulting from ingested anisotropic but homogeneous turbulence has been 
studied with existing fan noise models with particular emphasis on fan noise resulting from 
atmospheric turbulence in static and flight conditions. Based on the material presented in 
this study the following conclusions have been drawn and the following recommendations 
are made: 

Rotor Tone Noise Due to Atmospheric Turbulence 

• The difference in fan tone sound power at BPF due to ingested atmospheric turbulence 
between typical static test conditions and typical landing approach conditions is in 

the order of 30 dB. Based on this difference, it is concluded that fan noise due to 
ingested atmospheric turbulence, is negligible in landing approach conditions for the 
currently-used high-bypass ratio engines. 

• Fan noise due to atmospheric turbulence is expected to be below actual fan noise in 
flight conditions. It is therefore assumed, that it is not necessary to simulate in-flight 
turbulence in static tests. However, fan noise due to atmospheric turbulence in static 
tests should be reduced to levels below broadband fan noise in flight conditions. An 
effective way to achieve this is the reduction of the turbulence intensities in the static 
inflow field. The following guidelines are recommended for static test conditions: 

• usage of an inflow control device 

• low limits on wind velocities 

. small surface roughness scale in the test stand environment 
. large height of engine axis above ground 

. mean wind direction from the forward arc to eliminate reingestion 

Proper test stand and bellmouth inlet design and possibly inlet boundary layer suction are 
necessary to achieve this reduction. 

• Fan tone noise due to convected turbulence is dominated by the distortion elements 
within a small range of transverse scales. For typical turbulence energy distributions 
this range covers about one decade and is centered around a transverse scale in the 
order of 30% of the rotor blade spacing at the blade tip. For a given turbulence vari- 
ance, maximum sound power levels are achieved at transverse integral scales in the 
order of 25% of the rotor blade spacing at the blade tip. For a JT9D fan this is about 
4. 1 cm and for the JT 1 5D fan it is about 1 .3 cm. 
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• Although fan tone noise due to convected turbulence is dominated by the distortion 
elements within a small range of transverse scales, it is relatively insensitive to changes 
in the transverse integral scale. Over a large range of transverse integral scales (5 to 
200% of the fan blade spacing at the blade tip) the sound power due to convected 
turbulence varies only by 6 dB. This is due to the turbulent energy distribution, which 
is based on the von Karman spectra. For the mentioned region of transverse integral 
scales the range of distortion elements that dominate the fan noise generation contain 
always a significant fraction of the total turbulence energy. 

• All correlation coefficients characterizing the inflow turbulence field have an effect on 
the resulting fan noise. A vortex field with its uneven radial velocity distribution 
function causes lower sound power levels than the corresponding Gaussian eddy field 
with an even radial velocity distribution function. This is due to the fact that a 
significant amount of the turbulent energy in the vortex field contributes to higher 
order radial modes which do not propagate. 

• The difference in fan tone noise due to atmospheric turbulence between static and 
flight conditions is primarily due to the difference in the intensities of the turbulence 
involved in the fan tone noise generation. In static conditions the atmospheric turbu- 
lence is affected by a large inflow contraction, which results in: 

. large initial scale and therefore high initial turbulence energy level 
• large increase in the standard deviation of the transverse turbulence component 
. large streamwise integral scale at the fan face. 

In the flight conditions the turbulent energy at large wavenumbers decreases with 
altitude, due to the nearly constant turbulence variance and the increase of the integral 
scales with altitude. 

• Due to the mean velocity profile in the ambient air in static conditions, the highest 
intensities for the transverse turbulence component are expected to occur in the bot- 
tom of the inlet, where the turbulence is exposed to the largest contraction ratios. It is 
recommended that special attention be given the attenuation of the turbulence enter- 
ing the fan in the bottom segment of the inlet. 

• For the fan geometry representative for the currently used high bypass ratio engines, 
the fan noise due to atmospheric turbulence is due to the streamwise turbulence 
component in the flight condition and due to the transverse component in static tests. 

Effect of flow contraction on initially isotropic turbulence: 

• According to analytical models, the transverse turbulence component increases while 
the streamwise component decreases in a flow contraction. For small contraction ratios 
the analytical model predictions agree reasonably well with the experimental data. But 
at large contraction ratios, typical for static tests, the analytical models overpredict the 
attentuation of the streamwise turbulence component. 


133 



• Most of the changes in the spectral distribution of the transverse turbulence component 
occur in the first segment of the contraction, at contraction ratios smaller than three. 
The most significant change occurs in the one-dimensional spectrum for the transverse 
wavenumber component normal to the turbulence velocity component. For large 
contraction ratios the post-contraction turbulence field is dominated by the trans- 
verse turbulence component and can be represented with a unidirectional vortex field. 

Inlet angle of attack effects on fan noise; 

• The effect of inlet angle of attack on fan noise has been investigated with sound pressure 
measurements on the inlet walls, fan duct walls and fuselage surface over the full 
operational range of aircraft angle of attack. The measurements indicate that there is 
no obvious correlation between inlet angle of attack changes and changes in broadband 
noise as well as narrowband noise at BPF and its second and third harmonic. A circum- 
ferential variation of the buzzsaw sound pressure field in the inlet, that increases with 
angle of attack has been observed. The acquired data did not allow the evaluation of 
inlet angle of attack effects on fan noise radiation. The feasibility to evaluate such an 
effect in flight tests should be investigated prior to future flight tests. 

• For the particular engine used in the flight test the most significant changes in the 
mean flow field at the fan face, associated with changes in the inlet angle of attack, are 
probably due to a change in the rotor tip flow rather than changes in the inlet flow 
field. Most likely they result from the relative motion between the fan and the fan 
casing introduced by the asymmetric loads on the inlet nacelle. 

• The remaining fan tone noise at BPF in flight conditions is neither due to atmospheric 
turbulence nor due to inlet angle of attack related distortions. 
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APPENDIX 

ONE DIMENSIONAL POST-CONTRACTION SPECTRA 


Ribner and Tucker (Ref. 17) derived the three-dimensional post-contraction spectra for an 
axisymmetric contraction and initially isotropic turbulence 
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with G(k) = spectrum function of isotropic turbulence. 

Based on these equations and a given spectrum function, the one-dimensional spectra at the 
post-contraction location can be determined by integration over two of the wavenumber 
components. This integration is only possible, however, for simple spectrum functions. 

One of the simplest functions used to represent turbulence at large Reynolds numbers is 
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Based on this spectrum function Ribner and Tucker derived the longitudinal one-dimensional 
spectra at the post-contraction location 

F2®(kl)=F3®(kl) 

In the course of the present investigation three additional spectra have been derived to 
complete the matrix of one-dimensional post-contraction spectra 


F2®(k2) = F3'^(k3) 
F2®(k3) = Fs^lkj) 


Explicit equations have been derived for these spectra 
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Since the three-dimensional spectra have been formulated as a function of k rather than 
3C , the integration will be executed in the k coordinate system 



The integration is performed in polar coordinates 

k| = r cos 6 k3 = r sin 0 
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The integrations provide the following result 
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S = p+ 1 


(A. 15) 
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(A. 16) 
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This solution is restricted by the following condition 
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The one-dimensional spectra for the transverse turbulence component are based on the 
three-dimensional spectrum of equation (A. 2) and are obtained by integration over kj, and 
either k 2 or k 3 - Integration over k] and k 3 gives 
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s, q, A and e are defined by equation (A. 1 5) through (A.l 8 ). The validity of the solution 
is again limited to 




1 - € 


For large contraction ratios (e « 1 .) the spectrum can be approximated by 
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Integration of the three-dimensional spectrum defined in Equation (A. 2) over kj and k 2 
provides the second one-dimensional spectrum for the transverse component 
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(A.31) 


s, q, A and e are again defined by equations (A. 15) through (A. 18) and the validity of the 
solution is again limited to 




The approximate solution for large contraction ratios (e« 1 ) is 


The integral scales are determined based on the following relationship. 

' 1 ^ . -i Jfi Xj 

J ‘‘’‘j ^ 

_ J “ ■ 


= 7^ Fi® (Kj = o) FjB (kj = .0) 


The post-contraction integral length scales have been determined based on this equation. 
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F2® (k2 = -0) = 
F2® (k3 = O) = 
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= variance of precontraction turbulence 

= integral scale parallel to turbulence component in precontraction turbulence 
The ratios j and defined by equation (4.2.20) and (4.2.21). 
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